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Abstract 

Using the replica method, we analyze the mass dependence of the QCD3 partition function in a 
parameter range where the leading contribution is from the zero momentum Goldstone fields. Three 
complementary approaches are considered in this article. First, we derive exact relations between 
the QCD3 partition function and the QCD4 partition function continued to half-integer topological 
charge. The replica limit of these formulas results in exact relations between the corresponding 
microscopic spectral densities of QCD3 and QCD4. Replica calculations, which are exact for QCD4 
at half-integer topological charge, thus result in exact expressions for the microscopic spectral 
density of the QCD3 Dirac operator. Second, we derive Virasoro constraints for the QCD3 partition 
function. They uniquely determine the small-mass expansion of the partition function and the 
corresponding sum rules for inverse Dirac eigenvalues. Due to de Wit-'t Hooft poles, the replica 
limit only reproduces the small mass expansion of the resolvent up to a finite number of terms. 
Third, the large mass expansion of the resolvent is obtained from the replica limit of a loop expansion 
of the QCD3 partition function. Because of Duistermaat-Heckman localization exact results are 
obtained for the microscopic spectral density in this way. 
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1 Introduction 



QCD in 2+1 space-time dimensions (QCD3) provides an interesting alternative to spontaneous chiral 
symmetry breaking in 3+1 dimensions. Because of the absence of a 75 matrix in any odd number 
of space-time dimensions, ordinary chiral symmetry cannot be defined. Yet, for an even number of 
flavors 2Nf it has been known for some time that there exists a close analogue of chiral symmetry in 
2+1 dimensions (see, e.g., ref. |ffl]). In this case the global flavor symmetry may break spontaneously 
according to U(2Nf) — ► U(Nj) x U(Nj) |l], |2j. For an odd number of flavors it has been argued that 
the pattern of spontaneous symmetry breaking pattern is replaced by U (2Nj + l) — » U (Nf + 1) x U(Nf) 
Q. This picture of spontaneous flavor symmetry breaking has been supported by recent Monte Carlo 
simulations ||. 

Spontaneous flavor symmetry breaking such as discussed above gives rise to one Goldstone boson for 
each broken generator. In ref. [§] the low-energy effective partition function of the corresponding chiral 
Lagrangian was used to extract information on the spectrum of the Dirac operator of QCD3. It was also 
suggested M that Random Matrix Theory may be used to exactly compute the microscopic properties 
of the Dirac spectrum in much the same way as was first proposed for QCD in 3+1 dimensions (QCD4) 
0. The resulting spectral correlation functions for the ensemble relevant for QCD3 (a variant of 
the Unitary Ensemble, UE) are indeed universal in Random Matrix Theory |^, ||, M|. They can be 
expressed entirely in terms of (2+l)-dimensional finite- volume partition functions Q. Recently, the 
patterns of flavor symmetry breaking for theories with gauge group SU(2) or with fermions in the 
adjoint representation of SU(iV c ) in 2+1 dimensions have been shown to correspond to two other 
classes of the classical Random Matrix Theories [||, namely the Orthogonal Ensemble (OE) and the 
Symplectic Ensemble (SE), respectively, in precise analogy to the (3+l)-dimensional case pi. The 



microscopic spectral densities of these theories were calculated recently in |TT], 12 1 



To derive non-perturbative analytical results for the smallest eigenvalues of the Dirac operator it is 
important to prove that universal Random Matrix Theory results coincide with the low-energy limit of 
the field theory. For QCD4 this has been achieved for the spectral density and the two-point spectral 



correlation function by means of the so-called super symmetric method [TL3, 14, 15, M, 17, uM, 19, B 



In this method, one adds a set of fermionic and bosonic ghost quarks (or "valence quarks", in a 
terminology borrowed from lattice gauge theory [21, 22, ^3|) such that their determinants cancel for 
equal masses. In this way one obtains a generating functions for the resolvent (the partially quenched 
chiral condensate in the language of lattice gauge theory) or the higher order correlations functions. 
The spectral density then follows from the discontinuity of the resolvent in the complex mass plane of 
the ghost quarks |R 



An alternative method to eliminate the determinants of ghost quarks is the replica method [24]. 
It is based on adding either n fermionic or n bosonic valence quarks (also called replicas), performing 
an analytical continuation in n and taking the limit n — > at the end of the calculation. Although the 



replica method has been challenged [25], this method has received a considerable boost in the context 



of condensed matter physics during the past year [26, 27, 28, |29fl. Recently, the replica method has 



also been applied to the field theory corresponding to the low-energy limit of QCD4 |5Cj, 31], p^] . The 
advantage of the fermionic replica method is that we are familiar with the pattern of spontaneous 
chiral symmetry breaking, and it is straightforward to define the corresponding low-energy field the- 
ory. Although we have less experience with spontaneous symmetry breaking in the supersymmetric 
formulation, the consensus is that such internal supersymmetries are not spontaneously broken. Based 
on this assumption there should not be any problems in defining the low-energy limit of the super- 
symmetric partition function of this theory, which in turn is the generating function of the resolvent. 
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This has recently been confirmed by the explicit calculation of Szabo [ 33 1 , who has demonstrated how 
the microscopic spectral density and other spectral correlation functions of QCD3 can be obtained in 
this way. Nevertheless, precisely in the case of QCD3 it is advantageous to explore also the replica 
method. The reason is as follows. In general, one expects that the replica method can only provide 
us with asymptotic series for spectral correlation functions of the Dirac operator, and not with exact 
results. In QCD3, however, the large mass asymptotic expansion of the resolvent terminates (it is 



semiclassically exact |29] in the quenched case) and exact results will be obtained this way. 



An important ingredient of this paper is the discussion of general relations between the finite- 
volume partition functions for QCD in 2+1 and 3+1 dimensions. They will be introduced in section 
2 and the relations will be formulated and proved in section 3. These relations hold in all generality 
and are independent of any detailed derivations in any of the two theories. The implications of these 
relations are discussed in section 4, where we use the replica method to derive our main result: two 
completely general relations between the microscopic spectral densities of the two theories. Thus, once 
we know the spectral densities of QCD4, the spectral densities of QCD3 follow. In section 4.1 it is 
shown that results for the microscopic spectral density and microscopic spectral correlation functions 
obtained in this way are in exact agreement with previous results form Random Matrix Theory 
An exact calculation of the resolvent of the QCD3 partition function from its finite volume partition 



function is thus possible by using the recent results for the QCD4 finite volume partition function [31] 
(see section 4.2). In sections 5 and 6 we perform explicit calculations in QCD3 based on the replica 
method. First, in section 5, a small-mass series expansion of the resolvent of the QCD3 Dirac operator 
is generated by means of Virasoro constraints on the finite volume QCD3 partition function. As a 
by-product of this analysis we derive a series of spectral sum rules for the Dirac operator in QCD3. 
In section 6, we explore the (large-mass) saddle-point expansion of the resolvent, which turns out to 
truncate after a finite number of terms, and hence is exact. This truncation is consistent with the 
general relations to the QCD4 partition function derived in section 3. Our results are summarized in 
section 7. 



2 The QCD 3 Partition Function 



The Dirac operator of the QCD3 partition function is given by 

D = i(d k + iA k )a k , (2.1) 

where the a k are the Pauli matrices and the are the SU(N C ) gauge fields. The partition function 
for 2Nf flavors is defined by 

Z = (det(Dl 2Nf +iM)}, (2.2) 

where M is a Hermitian mass matrix and the average is over the gauge field configurations weighted by 
the Yang-Mills action. Without loss of generality, the mass matrix can always be chosen diagonal, and 
we will only consider this case. Generally, the fermion determinant in QCD3 is not positive definite, 



resulting in a Chern-Simons term due to the anomaly [34]. However, for quark masses occuring in 
pairs rrifc and —rrik, the fermion determinant is positive definite, and the theory is anomaly free. In 
that case the sign of the {tpk^k) is determined by the sign of the quark mass, 

(^) = -7^S. (2.3) 
\m k \ 
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The low-energy effective partition function is constructed from the requirements that its trans- 
formation properties should reflect those of QCD3, i.e. the partition function is invariant under 
U £ U(2Nf) transformations of the quark fields if, at the same time, the mass matrix is transformed 
according to 

M -> UMU- 1 . (2.4) 



The chiral condensate given by (2.3) is only invariant under an U(Nf) x U(Nf) subgroup of U(2Nf). 
The Goldstone manifold is thus given by U(2Nf)/U(Nt) x U(Nf). The matrices in this manifold 
transform according to 

U -» UiUU^ 1 , Ui € U{2N f ). (2.5) 

To lowest order in the mass, the mass term of the effective theory is given by 

C m = ±Re£TrM?7. (2.6) 

Because of the integration over the coset in the partition function, the overall sign of the mass term is 
irrelevant. For an even number of flavors in the limit where the Compton wavelength of the pseudo- 
Goldstone bosons is much larger than the size of the box, the partition function factorizes into a zero 
momentum part and a non-zero momentum part. The first one, which is also known as the QCD3 
finite volume partition function, is thus given by 

ZqSdAM) = I dUexp[VZTr(MUT 5 U^] , (2.7) 

w 3 JU(2Nf) 

with A4=diag(mi, . . . ,mjvy, —fnii ■ ■ ■ , — ^iv* )Q, r5=diag(ljv / ,-lj\r f ) and V the volume of the three- 
dimensional space-time. We have used an explicit representation of the coset, U — ► UT^W, and 
extended the integration to the full U(2Nf) group. This gives only rise to an overall volume factor. 



The partition function (2.7) is only a function of the rescaled masses 

IH = rmVZ, (2.8) 

which will be kept fixed as V — > 00. We shall here consider QCD3 at fixed ultraviolet cut-off A. In the 
infrared the theory is regularized by the three- volume V, which is only eventually taken to infinity. 
If not only the volume is finite but also space-time is discretized such as in lattice QCD, the total 
number of Dirac eigenvalues is finite and will be denoted by J\f. 

For an odd number of flavors one of the quarks is unpaired, and the fermion determinant may 
change sign under large gauge transformations. In this case the sign of the chiral condensate is not 
given by the sign of the quark mass. There are two inequivalent possibilities for the vacuum state, 
determined by the sign of the fermion determinant. They are not connected by unitary transformations 
of the coset elements. In a regime for which m ~ 1/VS, both states have to be included in the low 
energy limit of the QCD partition function. The QCD3 partition function satisfies the parity relation 

ZQCD 3 (-m) = (-l)^-Z Q c D3 (m). (2.9) 
For an odd number of flavors the low-energy effective partition function is thus given by 

zg^ttM) = / dUcosh[VXTr(MUf 5 tf)} , (2.10) 

yuu 3 Ju{2N f +l) 

1 Recall that in 2+1 dimensions there is no definite sign for the fermionic mass term, a consequence of the fact that 
in any odd number of space-time dimensions the two sets of 7- matrices {7^} and {— 7^} are inequivalent irreducible 
representations of the Clifford algebra. 



4 



for even number of eigenvalues M, and 

Z { ^ 1] ({^}) = [ dU S mh[VXTr(MUf 5 U^)} , (2.11) 

^ u 3 Ju(2N f +l) 

for odd number of eigenvalues. The second partition function (2.11) is not positive definite, and is 



unsuitable as a physical partition function. In particular, it vanishes when the unpaired fermion mass 
vanishes 0. Nevertheless, this second partition function can appear at intermediate stages Q. 

These partition functions are to be compared with the QCD4 finite volume partition function given 
by i 

2& N,) (.M) = I dU (det U) v exp \]-VT,Tt(MU^ + UM*)\ , (2.12) 
JU(Nf) L2 

where v is the topological charge, and ,M=diag(mi, . . . , rnjsr,). In both cases £ stands for the infinite- 
volume quark-antiquark condensate. As we will see in the next section, there are deep relations 
between these two types of group integrals. Although the structural nature of these relations suggest 
that they can be obtained from general group theoretical arguments, we have only been able to prove 
them at the technical level. 

3 The QCD3-QCD4 Connection 



In this section we shall prove a set of surprising relations between group integrals of Itzykson-Zuber 
kind (2_ 1 7), ( |2. 10 ) and ( 2.11| ), and group integrals of "external field" type (|2 . 1 2|) . Some of these relations 



were already pointed out in ref. 35], based on the relationship to Random Matrix Theory, but we 



shall here prove them directly from the finite volume partition functions. 



Both of the group integrals in eqs. ( |2.7D and ( p. 12 ) can be expressed in closed form. The first is 



simply an Itzykson-Zuber integral |36| with a particular iVj-degeneracy in one of the matrices (T5) 
and the result is || 

ie J A({ ft }) A({- ft () \ 

agK«*» = ( -i)W>/' U«-*>> ) _ 



where the NfxNf matrix A({/Xj}) is defined by 

A({ tt }) il = (ftfVi , j,l = l,...,N f . (3.2) 
and the Vandermonde determinant is taken over all masses (including the sign- mirrored ones), 

2N f 

A(M) = Him-m) ■ (3.3) 

i>3 

Recently, it was shown that the above partition function can be expressed much more compactly by 
means of just a single Nf x Nf determinant ||33|| . However, precisely for our present purposes the 
above form is actually more convenient. 



2 Note that we have chosen a different prefactor than in ref. || in order to make the partition function positive 
definite. For averages the difference is of course irrelevant. 
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For an odd number of flavors the integral ( |2,10| ) and ( 2.11| ) can also be evaluated as Itzykson-Zuber 
integrals. The result is 



r (2JV>+l), 
•'QCD± 



2 N t 1 r 



detD( At ,{// i })±(-l) JV/ detD(-//,{-/i i }) , (3.4) 



A(M) 2 

where the (2Nf + 1) x (2Nf + 1) matrix D is defined as 



2"/ detDO*, M) ^ detf A (MMk /+ ix;v /+ i A(-/x , {-^}) N/+1xN/ \ (3 _ g) 



Here, the matrix A is given by eq. ( |3.2|) , with the corresponding quadratic or rectangular size explicitly 
indicated. 



For integer v the integral ( 2.12 ) can be expressed as [37|, 

d N, \{lH}) = detB({^})/A({^ 2 }), 
where the matrix B in eq. ( |3.6|) is given by 

B ({M)j« = Mji -1 ^ -1 ^ (Mi) , j,l = l,-..,N f 



(3.6) 



(3.7) 



with Iu = (d/dfj,) 1 The denominator is given by the Vandermonde determinant of 

squared masses, 

N f 

A(K 2 }) = II(Mi-M?) = detf^- 1 ] . (3.8) 

In both cases the normalization convention of the integrals has been chosen for later convenience. 

For non-integer values of v we define the QCD4 partition function according to the analytical 
continuation in v of the Bessel functions iBTf, 



e zco S e e iue de _^a 



7T JO 



dse 



-z cosh s—vs 



(3.9) 



and not by analytical continuation in v of the unitary matrix integral ( |2.12 ). 

The integral ( 2.12j ) satisfies an important relation known as flavor-topology duality p8[. A proof 
that greatly simplifies earlier proofs [38] is obtained by realizing that the matrix B in ( [3.6] ) can be 
replaced by 

B({M)ii -> M} -1 -W-i(Mi) > j,l = l,...,N f , (3.10) 

If we take the limit /U — ► and use the fact that I v (n) ~ (n/2) u /T{v + 1) we find the the flavor-topology 
relation 

. 1 : (3.11) 

,=1 ^ 

The general relation obtained by iterating ( |3.11| ) with repect to Nf is valid for arbitrary v provided 

A N f) ; 



AT* 



that -2^ is analytically continued in 1/ according to Q3.9| ). 

After these preliminary definitions, we are now ready to state and prove three theorems. 
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Theorem I - Even number of flavors in QCD3: 

Let Zqcd {{Pi}) an d Zv* ({/^i}) be as defined in eqs. ( \2. TQ and with the normalization 



conventions 1\) - (\3.{\). Then the following identity holds: 



Proof: We start by simplifying the right hand side, using that for v = ±1/2 we have 

4=-i/2<» = \ — cosh^u) and ^=+i/ 2 (^) = J — sinh(//) . (3.13) 

y ' y 7T/i 

One can easily convince oneself that in eq. (^^) only terms with all derivatives acting on cosh (//.,•) 
and sinh(/xj), respectively, contribute to the partition function. All other terms containing derivatives 
of the square roots can be eliminated by adding multiples of columns from the right, which leaves the 
determinant invariant. We can thus pull out the common factors of square roots and arrive at 

N, 



Z { Z% 2 {{^}) = U\— detC({M)/A({M?}), C^M^M^cosh^ 1 ) (/,,■) 
j=i V ^ 



f / r\ 

= II \ — detS({M)/A(K 2 }) , Sa^^^sinh^ 1 )^). (3.14) 
j=i V 

Here and below cos or sinh^(/u) denotes the p-th derivative of cosh(mu) or sinh(mu). Turning 

next to the QCD3 side, the Vandermonde determinant can be written as, 

N f 

A(M) = [] 2m* A 2 ({/, 2 }) . (3.15) 

k=l 



In order to bring the (2Nf) x (2Nf) determinant in eq. ( |3.l[) into the desired form we iteratively add 
and subtract columns to obtain cosh and sinh. Starting with the first column we add the Nf-ih column 
to obtain (2 cosh(//i), . . . , 2 cosh(fj,N f ), 2 cosh(^i), . . . , 2 cosh^A^ )). Then, subtracting 1/2 of this new 
first column from the Nf-ih column, we obtain (— sinh(/ii), . . . , — sinh^A^), sinh(^i), . . . , sinh^TVf )) 
as a new column vector at position Nf. The common factor of 2 can be taken out of the first 
column. We proceed in the same way with the second and (Nf + 2)-nd column, where cosh and 
sinh get interchanged in the resulting columns compared to the previous step. After completing this 
alternating procedure we arrive at the following block structure for the QCD3 side: 

z^({ u \) - ( i)Nf(N f +m 2Nf dct I C ({M) - S ({MA 

^qcd 3 UM) - ( i) A(M) det [c({^}) aOta})) 

A(.M) V C({/ii}) S({^}) J 

where in the last step we have added the lower blocks to the upper ones, by a subsequent addition of 
rows. We can now perform a Laplace expansion into products of Nf x Nf blocks, which due to the the 



zero block, 0, only yields one single product of two Nf x Nf determinants. Together with eq. ( 3.15 ) 
we obtain the final result: 



2qS(W) = — m detC({M)detS({M) • (3-17) 

A 2 ({^ 2 })n fc £i^ 
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By comparing eqs. ( |3.17 ) and ( |3.14 ) we can read off the prefactor in the theorem eq. ( 0.12 ). 



Theorem II - Odd number of flavors in QCD 3 : 

Let in addition Z^^^lp, {fa}) be as defined in eqs. (3.4) and ( \3.5] ). Then the following two 
identities hold: 

^\fa{fa}) = ^/f ZI%\%(» M ) zi% 2 ( M ) . (3.18) 



Proof: Since the right hand side already follows from eq. ( 3.14 ) we start with the QCD3 side, 
treating both cases, ±, simultaneously. The Vandermonde determinant with one additional mass can 
be treated similarly to relation ( 3,15| ) and we obtain 

N f N f 

ACM) = (-l)^^ +1 )/ 2 n(^-^)Il 2 W A(K 2 }) 2 • (3.19) 

j=l 1=1 

Next, we proceed by adding and subtracting columns as in the proof of the previous theorem. We 
arrive at the same determinant with blocks of matrices C and S, except that the (Nf + l)-th columns 
remains unchanged. Here, we are dealing with an odd number of flavors and the determinant in eq. 

(M) 



2 N t detT>(fa{fa}) 



2 N t det 



is of size (2Nf + 1) x (2Nf + 1), with leaves one column unpaired. We thus have for eq. ( |3.5| 

/ pi- 1 cosh^Gu) n N f& -fjP^smk^Xn) 

C({fa}) tf'eP S({fa}) 

V C({^}) (-m) N fe-" S({fa}) 



(3.20) 



When adding the second determinant of eq. (|3.4| ) with negative arguments we can use the fact that 
the matrix C is even whereas S is odd. After taking out all the minus signs from the matrix S, which 
gives (— 1) /, we can add the two determinants defD(fa {fa}) and det D(— fa {— fa}), as they now 
differ only by the (Nf + l)-th column. The resulting determinant reads 

(3.21) 



detD(MM) ± (-l) JV/ detD(-/i,{-/i i }) 

/ /^cosh^ 1 ^) 



det 



y-isinh^'-^Gu) 



fi N f(e"±(-l) N fe-") 
C({fa}) /i^(e«±(-l)^e-«; 
V C({fa}) ±^ f (e^±(-l) N fe~^ 

We now have to treat each of the two cases db separately. For the "+" case the (Nf + l)-th column 
reads (2fa N fcosh ( - N f\fa), 2fa f cosh^ N f\fa), 2fa f 'cosh^^i), . . .) which can be absorbed into the 



-s({M) 
s({M) 



matrices C. The "-" case yields (2fa N fsmh ( - N f ) (n), 2fa f smh^ N f\fa), -2/4 'sinh^Vl^ ■■■), and 
we can absorb it into the matrices S after changing the sign of the vector and permuting it to the 
very right. We will then add or subtract blocks again to produce a Ojv> x jv> block on the right or left, 
respectively. In more detail we have for the 



det DGu, { W }) + (-1)^ det D(-M-M) 



"+" case 



2 det 



(3.22) 



2 det 



^cosh&'-^t/i) 
C({^i})jv>xJV/+l 
C({fli})N f xN f +l 

ji*- 1 cosh&'-^V) 

2C({/ij})Ar /X 7v / +l 
C ({fa})N f xN f +l 



S({l J 'i})NfxNf 
S({fa})N f xN f 

-/^sinh^V) 

ONfXNf 



s({M) 



NfXN 



f 
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Because of the zero block ON f xN fl a Laplace expansion into products of (Nf + 1) x (Nf + 1) blocks 
and Nf x Nf blocks contains at most Nf + 1 nonzero terms. However, all the determinants for which 
two rows of the matrices C occur twice, also vanish. Thus we are left with 

det D(/i, {fj,i}) + (-1)*/ det D(-/x, {-//*}) = 2 JV >+ 1 det C(/i, {/*}) det S({M) , (3.23) 

where the size of the determinants on the right hand side is given by the number of arguments. 



Along the same lines we obtain for the "-"-case, 



det D(/x, {m}) - (-l) N f det D(-//, {-^}) 



(-!)*/+! det I 



(-1)^+1 det I 



/^cosh^ 1 ^) 

C({fLi})N f xN f 
C({Hi})N f xN f 

ONfXN f 



sm 



h (i- 



(3.24) 



S({Hi})N f xN f +l 
^({^i})N f xNf+l 



-/'■ 



C({/Xi}) 



-2S({/ii})jVj. xiV/H-1 
S({Mi})jV)xiV)+ly 



2 7V /+ 1 detC({^})detS(^,{^}) 



In the last step the minus signs coming from permuting the (Nf + l)-th column of eq. ( 3.21 ) get 
absorbed into the sign of the only non-vanishing contribution from the Laplace expansion. Taking 
eqs. ( 3.23 ) and ( 3.24 ) together with eq. ( |3.19| ) we finally obtain 



„(2N f +l), , ^ 
Z QCD± (MM) 



(_l)iv/ 2 W / 



A({^}) 2 n£i(^ 2 -^)w 

We can now compare to the right hand side of eq. ( 3.1S| ) with the help of eq. ( |3.14 ). Looking at 

x| — /x 2 ) of the Vandermonde determinant of zj) Nf+1 1 
factor (—1)^', which adds up to the correct prefactor of eq. ( 3.1£| ) 



det C(n, {fii}) det S({^}) for + 
det S(/i, {fii}) det C({/^}) for - 



(3.25) 



the factor rifcG^fc ~ ^) °f the Vandermonde determinant of Z^ s ^ x (/i, {ni}) we obtain an additional 



The first two theorems stated so far relate the QCD3 partition function with an odd or even 
number of flavors to a single product of two QCD4 partition functions. In the following we wish to 
consider the QCD3 partition function with additional flavors which do not come in pairs of opposite 
sign. However, such partition functions do not satisfy relations as simple as Theorem I and II. Instead, 
they are given by linear combinations of products of QCD4 partition functions and do not directly 
give relations between the spectral density of QCD3 and QCD4. For this reason we will only treat the 
simplest case, that of the even flavor QCD3 partition function with two additional unpaired masses. 



Theorem III - Additional unpaired flavors in QCD3: 
Let in addition Zqqv ({jM,Ci w ) ^ e defined as 



det / A({M,C) A({-fn},-0 



A(M) 



(3.26) 



where the (Nf + 1) x (Nf + 1) matrix A({/ij}, Q is defined as in eq. (j3.$\ ). The mass matrix is given 
by A4 = diag(/ii, . . . , fJ>N f , C; ~~ Mi> • • • j ~l i Nti U} )- Then the two following relations hold: 

(3.27) 
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PROOF: In order to prove these two statements we will first rephrase them in terms of linear combi- 
nations of QCD3 partition functions^]: 

r (2iV/+2). 



and 



-QCD 3 G>i},C,^) + Z QCD 3 = 



(3.28) 



(3.29) 



Clearly, when taking the sum or difference of these two equations we find back the theorem eq. ( 3.27 ). 
In the following we will reduce eqs. (3.28) to identities between determinants which will be then 
proven separately in Appendix |A[ 

First of all, we note that, as in the proofs of the previous theorems on the QCD3 side, we can again 
reorder the determinant in terms of hyperbolic functions. After permuting the additional arguments 
C and ui we obtain 



(2JV /+2 ) a , . , )(Nf +1)(JV +2 )/2ojV f +1 det D({^}, C, hQ 



(3.30) 



where we define the (2Nf + 2) x (2Nf + 2) determinant D for an even number of flavors, 

/ C({/ii}) -8(0**}) \ 

,, n/ , u x , . C ({M) S({//i}) 

detD({^},C,w) = det C i-i cosh (i-i) (c) -(J-isinh^- 1 )^) ' ( } 

V a^cosh^H -(J^wH^-^iu) J 
In the next step, we write out the Vandermonde determinants explicitly as in eq. ( [3.19D , 

Nf N f 

a({^},{-m,±c,«) = (-i)^ ( ^ +i)/2 (^TC)n 2 ^ A (^ 2 }) 2 n( w2 -^)(c 2 -^ 2 ) (3-32) 



3=1 



for the QCD3 side, and 



Nf 



A(K 2 },C 2 ,- 2 )A(K 2 }) = (o; 2 -C 2 )A(K 2 }) 2 n(^-M|)(C 2 -^) 

i=i 



(3.33) 



for the QCD4 side. Putting everything together, we obtain from eq. (3.28) 



(u; + C)detD({ w },C,u/) + (w-C)detD({^},-C,w) = (-1) JV /+ 1 2 JV / +1 det C({M, C, det S({^}) 

(3.34) 

and from eq. ( |3,29| ) 

(^ + C)detD({M,C,w) - (w-OdetD^h-Cw) = (-l) N f2 N f +1 det S({ W }, C, w) det C({M) • 

(3.35) 

These two relations between determinants are derived in Lemma 1 ( |A.1| ) and Lemma 2 ( |A.1[ ) of 
Appendix A, which completes our proof of Theorem III. 



3 Note that there is nothing new in the definition ( |3.26| ) as compared with the original definition ( |3.l[ ). In order to 
avoid confusion, we have indicated explicitly where the unpaired masses should enter in the individual blocks. 
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4 Applying the Replica Method: General Results 



In this section we will use the replica method to derive general relations between the microscopic 
spectral densities of QCD3 and QCD4. 

Let us first consider QCD4 with Nf fermions in the fundamental representation. We add n valence 
fermions with degenerate (rescaled) mass \x v to this theory. In such a theory, the resolvent, or partially 
quenched chiral condensate, is defined as 

v = fim-^— 0«,{M) • ( 4 - 36 ) 

L, rc— >0 n OfJL v 

The partially quenched chiral condensate in QCD3 is defined similarly (recalling that here the number 
of flavors is 2n). The microscopic spectral density of the Dirac operator is given by the discontinuity 
of the resolvent across the imaginary axis [13], 

Ps\C;M) = 2^ Disc s "(/^>{w})L=ic 

= i-[S,« + e;{M)-^(iC-e;{M)] (4-37) 

The analysis of Dirac operator spectra by means of a generating function of the resolvent is not 
restricted to just the spectral density. The fe-point spectral correlation functions can be computed 
from partially quenched chiral susceptibilities xi^vu ■ ■ ■ , ^v k ,{fJ-i})- These quantities can be derived 
from the partition function by extending it with kn valence fermions in k mass-degenerate sets. Then 
the n — > limit is again taken after appropriate differentiations, 

X^(n Vl ,...,n Vk ,{iJ,i}) = lim-^-^ In zi Nf+kn) (/x^,..., /J, Vk , {m}) . (4.38) 

In a spectral representation, 

(k) { r n [°° P^coim.(Cl>---)Cfe;{Mi}) , . 

xl>(Hv 1 ,---,Vv k ,{Vi} = / dd-'-dCk—— — — -, 4.39 

J -00 «1 + il vx ) • • • «fe + Vv k ) 

and we can again invert this relation to get the spectral correlation function as the discontinuity across 
the imaginary mass axis, 



(Ci,...,a;{M) = 77^ Disc xl fc) (/v,---,^ fc ,{M) _., • ( 4 - 4 o) 



0) 

^5, conn, vs 11 • • • > SKj in J ) (27r) fc 



For instance, for the 2-point function the explicit relation reads after using the symmetry of the 
spectrum around the origin 

P^Lnn.(Cl>C2;{Mi}) = ^b<i 2) «i + e,i( 2 + e;{fii}) +x { „ ) (i(i -e,*Ca 

-X { u ] «l " e> <2 + e; M) - x { u ] «i + e, 1(2 - e; {/*})] • (4.41) 

A similar formula gives us /O^onn (Cli ■ ■ ■ 1 Ck] {^i}) for arbitrary k. 
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4.1 The Microscopic Spectral Density 



Let us first apply the relation ( [4.37 ) to the partition function identity ( 3.12 ), adding to the 2Nf 
physical fermions n valence fermions of mass fi v and n valence fermions of mass — fj, v . The number 
n is eventually taken to zero. We then obtain a relation between the microscopic spectral density of 
the Dirac operator in 2+1 dimensions for 2Nf fermions pairwise grouped with opposite mass and the 
corresponding microscopic spectral densities in QCD4 evaluated at half-integer topological charge, 



/4Si(C;{M) 



1 

2 

1 r 



^(C;{a*})+pSSX(-C;{a*}) 



(2JV/) , 



pSC +1/2) (C;{^»+^ 



(N f) u=-l/2) 



(C;{M) 



(4.42) 



In the first equality we have used that the average spectral density is an even function of (. Note that 
the /z-independent factors in eqs. ( gup and ([yJD are immaterial. The same relation can be derived 
from the Random Matrix Theory representation of the finite volume partition function (see Appendix 
B). 

For 2Nf + 1 flavors we only consider the case Z^^+ for which the spectral density is positive 



definite. From the replica limit of the resolvent applied to ( |3.18| ) we find 



^Q 2 sr (c;{M,/i) = 2 



{N f ,v= 
^QCD 4 



- 1/2 \C;M,ri + ^X +1/2 ' 



«;{»}) 



where we again have used that the average spectral density is an even function of Q. For fi 
relation can be simplified by applying flavor-topology duality [38|, 



N 



1 



7T 



z 



(Nf) 
v=+l/2 



(4.43) 
this 

(4.44) 



Here, the prefactor gets modified compared to eq. (|3.18| ) for the following reason. If we define all 
partition functions normalized according to eqs. ( |3,6j ) and (3/7) the flavor-topology shift for general v 
(see eq. ( |3 . 1 1| ) ) yields the following factor: 



lim 



(JV/+1) 



_1/2(M>{^}) 



1 t(JV>) 



II 7, . Z u 1 2 

3=1 ^ 



Eq. Q4.44D leads to 



Pqcd 4 (C; M) 



(4.45) 



(4.46) 



a relation first discovered by Christiansen Q in the context of Random Matrix Theory. Here we see 
that this relation follows from the effective field theory partition functions alone, after using the replica 
method. Note how the "topological" u/fi-teim (for v = 1/2) nicely has cancelled out on the right 
hand side. 

The relations flO|) and flOg) are our two main results. They show that the microscopic spectral 
densities of QCD3 follow directly from those of QCD4 by analytic continuation in v. 

It is straightforward to check the relations ( 4.42j ) and ( 4.43 ) against the well-known expressions for 
the massless microscopic spectral densities for the two theories, as for example computed in Random 
Matrix Theory M, |, 
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pgSl(C;te = o}) = C j 



JN f -l/2{() 2 ~ ^+1/2(0^-3/2(0 



pSS" 1} (C;{/ii = o},o) 



pSSg?(C;Oii = o}) = I 



+ ^a^+iMO 2 - ^+3/2(0^-1/2(0 

<^Af / +l/2(C) 2 - "^+3/2(0^-1/2(0 
JN f +u{() 2 - JN f +u+l(C)JNf+u-l(C) 



(4.47) 



The relations derived above are indeed seen to be satisfied. 



It is also instructive to write down the asymptotic large mass expansion of the resolvent of the 
QCD3 partition function from the known asymptotic expansion of the resolvent for the QCD4 partition 
function. We only consider the case of an even number of massless flavors. The starting point is the 
relation (|3.12| ) for 2Nf massless flavors and 2n valence flavors with masses in pairs ±/i„. Using flavor- 
topology duality [38], the QCD4 partition function can be rewritten using 



lim \ \ 



r(* + *)(- 



/ 2 y+fc-il {Nf+n) Z^JM) 

[-) Z» ({/V,/W) = n{v+Nf) 



(4.48) 



This relation is readily derived from the explicit form of the finite volume partition functions 
From (ETI^ ) 

we then obtain 



T N f +n 



(2N f +2n) ( , _ 7T 1 (n) ( n ) 

/it) 



(4.49) 



After taking the replica limit we find the relation 



J QCD 3 if^v) 



r N f =0 ( 
Zu v=N s -l/2 



+ 



f+1/2^) 



El, 



(4.50) 



The replica limit of the quenched QCD4 chiral condensate is given by (for — |-7r < arg fi v < ^) |3l| 



E v (fi v ) _ . i(-l) u e- 2 ^ + {Au 2 -l){l-i(-l) u e- 2 ^) i(-l) u e- 2 ^(4u 2 - \){Av 2 -9) 



8^ 



64/z^ 



(4z/ 2 - l)(4i/ 2 - 9) [i(-l) u e- 2 ^ (19 - 4z/ 2 ) - 6] 
+ 3 x 256/4 +••• • 

We thus find the following large mass expansion for the resolvent of QCD3 



(4.51) 



v 27V / t \ 

■^QCDoUM 



Nf N 



f 



N f (N 2 



1) 



+ e 



("I) 



V f 



4/4 

Nf N f (N 2 -l) N f (N 2 -l)(N 2 -3) 
2 M 2 2 M 3 4 /4 



+ 



(4.52) 



For Ay = 1 and Nf = 2 all higher-order terms vanish. A replica calculation thus provides us with the 
exact non-perturbative result. 
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4.2 Higher Order Correlation Functions 



Let us now apply the relations between partition functions in QCD3 and QCD4 to higher order cor- 
relation functions. Because the valence masses occur in pairs ±/ij. in the QCD3 partition function, 
differentiation with respect to the valence quark masses does not exactly result in the Appoint spec- 
tral correlation functions. Rather we obtain nontrivial relations between correlation functions with 
arguments that differ by a minus sign. As an example, let us look in more detail at the two-point 
correlation function. We find 

PQCD 3 ,conn.(Cl 5 C2;{M) + PQCDgXonn. (&! "C25 {&}) 
+ PQCD 3 ,conn.( _ Cl,C2;{Mi}) + />QCD 3 ,comi. ( _ Cl, ~C?, {fii}) 

= p& 1 . /2) (Ci,C 2 ;{w}) + pt = om 1 . /2) (Ci,C 2 ;{w}) - (4.53) 

With connected two-point correlators given by the square of the kernel, 

^SLconn.(Cl,C 2 ;{/xa) = - (fo}, Cl , ft' " 



pgLo.(Ci,6;{M) = -[K^\M,<l(i)f , (4.54) 

the relation ( [4.53| ) is in complete agreement with result from Random Matrix Theory obtained in 
Appendix B (see eq. ( B.lOj) ). Similar relations can be derived for higher order correlation functions. 



We stress that the above relations are totally general, independent of any detailed evaluation of the 
n — > limit of the partition functions involved. The replica method thus provides us with truly non- 
perturbative relations which go beyond asymptotic expansions. Actually, in this case the asymptotic 
expansions terminate, as we shall see below. 

Finally, a note of caution for using the replica method. In our proofs we have not excluded the 
possible appearance of terms that vanish for positive integer values of n but are finite for n —* (for 
example, terms of the the form sinn-7r/n). Because of the agreement with the exact results we know 
that such terms are absent. It is an interesting problem to understand why this is the case for an 
analytical continuation in the number of flavors whereas such terms appear, as we have seen in eq. 
(|3.9|), when we perform an analytical continuation in the topological charge. 

5 Virasoro Constraints and Spectral Sum Rules 



In the next section we will turn to explicit calculations based on the replica method. So far there are 
no known techniques that have succeeded in taking the replica limit n — * beyond series expansions. 
The most obvious expansion that comes to mind is a Taylor series in the mass. Such an expansion 
was recently considered in the context of QCD4 where it was found that the so-called de Wit— 
't Hooft poles prohibit the evaluation of the partially quenched chiral condensate Yi u {fi v , {fJ>i}) beyond 
a given order, depending on the value of Nj and v. In QCD4, these de Wit-'t Hooft poles occur at 
precisely the order for which the logarithmic terms appear in the expansion [30, As a particular 



consequence, the small-mass expansion cannot be used to derive spectral correlation functions. There 
is simply no discontinuity across the imaginary mass axis up to the order at which the replica method 
predicts the partially quenched chiral susceptibilities^ 

4 Resummations of the series may hold a clue to how this obstacle can be overcome; see section 4.2 of ref. [pl[ . 
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As a preparation for the next section we derive in this section the small-mass expansions of the 
partition functions ( |2,7D and ( |2.10| ), Taking advantage of the results obtained we also derive the 
corresponding spectral sum rules, thus generalizing previous results. The small-mass expansions of 
( p.7| ) and ( 2.10| ) can be obtained by several different methods such as the character expansion of p6|]p|. 



We can also use the known results for the effective partition function of QCD4 as a generator for 
unitary integrals (as the effective QCD4 partition function is simply an external field integral over 
the group U(iVj)). However, the most elegant and efficient way to push the small mass expansion 
to arbitrarily high orders is probably through the use of Virasoro constraints. The appearance of 



Virasoro constraints in unitary integrals was first noticed in [41] for the case of large masses. It was 
later generalized for the case of small masses in [42]. In both works, only the one-link integral ( |2.12| ) at 
zero topological charge {y = 0) was studied, although in the small-mass expansion the flavor-topology 
duality readily gives the analytical continuation to v ^ 0. 

For the derivation of Virasoro constraints for QCD3 it is convenient to start with the following 
more general partition function: 



2f'(M) = J dUexp [TrMUIqU^ , (5.1 



where M is an arbitrary N x N matrix and I q is such that I? =1nxN and Trig = q . In particular, 
the even-flavor partition function of eq. ( |2.7| ) is obviously of this kind (with q = and M = VT,Ai), 
and it turns out that in the small-mass expansion other detailed properties of I q are irrelevant. Thus, 
in the small- mass expansion the partition functions ( |5.1|) with q = 0, ±1 are the building blocks for 
the QCD3 partition functions given by Q2.7p , ( |2.10D and (2.11). It was noted in ref. []30| that the 



finite-volume effective field theory partition functions satisfy the same type of differential equation, 
whether in 3 or 4 dimensions. Indeed, from the properties of I q we find 

/ d 2 



dM ca dM ab 



5cb 


)4 N) 




dZ™ 




dM ab 



= (5.2) 

= qZW (5-3) 
Due to the unitary invariance of the measure it is clear that the partition function, (|5.1f), depends 



only on the traces of the matrix M. Thus, following [42], we can change variables from the matrix 
elements M ab to 

t k = jTrM k , k > 1. (5.4) 



Clearly, for finite N not all tk are independent. Using the chain rule, the differential equations ( |5.2| ) 
and (|5.3D can be rewritten as 

J2 M£T 2 [£, - <M Zf ) = 0, (5.5) 
(A-g)zf) = 0, (5.6) 

where the operators C s are defined by 

s-l 

£ S = J2 d o d s-j + E ih d j+s + N8 S , 8>1. (5.7) 

3=1 3>1 



5 We have explicitly verified that the more general expansion we shall provide below agrees up to fourth order with 
the character expansion of |3q] for the case Tr(M) = 0. 
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We use the notation dj = d/dtj, and our convention is that for s = 1 the first term in £ s is absent. 
The £ a satisfy the Virasoro algebra^] 

[£ r ,£ s ] = (r - s)£ r+s . (5.8) 



A sufficient condition for the solution of the differential equations ( |5.2[ ) and ( |5.3| ) is thus given by 
so-called Virasoro constraints: 

£ s zf ) = (S s>2 + q6 aA )zW , 8>1. (5.9) 

They differ from the Virasoro constraints for QCD4 by the presence of the term 5 St2 and the parameter 
q which takes the value q = 1 in the QCD4 case. 

We normalize the partition function by Zq N \tk = 0) = 1. The Virasoro constraints can then 
be solved iteratively order by order in the masses by letting i& = after acting with the differential 
operators. It is instructive to consider the solution of the lowest few orders. Using the Ansatz, 

= 1 + ah + htl + b 2 t 2 + 0(/U 3 ) , (5.10) 



we deduce from £\Zq = qZg N ^ at tk = 0, 



Na = q. (5.11) 



At order we can use di£\Zq N ^ = qd\Z\ N ^ and £ 2 Zq N ^ = Zq N ^ at tp. = to obtain the equations, 

q 2 

2Nb 1 + b 2 = qa = ^, 

2b x + Nb 2 = 1. (5.12) 

The solution is given by b 2 = (N 2 - q 2 )/(N(N 2 - 1)) and 61 = (q 2 - 1)/(2(A 2 - 1)). Proceeding like 
this, we have derived the small-mass expansion of Zq N ^ up to order [i®, and the expansion can easily 
be pushed higher. We quote here the 0(/x 4 ) result, which is given by 

(m q (A 2 — q 2 )TiM 2 (q 2 - l)(TrM) 2 2q(N 2 - q 2 ) , 
Z^ N i = 1 + -f-TrM + V — — + W — ^ — ^ ttttt ;rTrM 3 
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A 2N{N 2 - 1) 2(N 2 - 1) 3iV(iV 2 - 1)(A 2 - 4) 



+ ^ TrMTrM 2 + 

+ 2(iV 2 -l)(A 2 -4) + 6iV 



N 2 {q 2 -3) + 4-2g 2 
(A 2 -l)(iV 2 -4) 



(TrM) 



3 



(iV 2 — q 2 )(5q 2 + 4 — A 2 ) ^ . (A 2 — g 2 )(3A 2 + 3g 2 — 2A 2 g 2 — 12) o 

+ — — — ' — TrM + ' — -TrMTrM 3 

^ 4A(A 2 -l)(A 2 -4)(7V 2 -9) 3iV 2 (7V 2 -l)(7V 2 -4)(7V 2 -9) 

(A 2 - g 2 )(24-10A 2 + iV 4 -6 g 2 -AV) (N 2 -g 2 )(q 2 -l) ^ 

8A 2 (A 2 -l)(7V 2 -4)(iV 2 -9) 1 ' 4A(A 2 -l)(A 2 -9) 1 ' 

4- g 2 [^V-6)-8A 2 ( g 2 -5)+6( g 2 -4)]+3iV 2 (iV 2 -4) 

+ 24iV 2 (A 2 -l)(A 2 -4)(A 2 -9) (lrMj + U(M 

In Appendix C we give the expansion up to order /i 6 . At order /u 5 we have found simple poles^] at 
A 2 = 0, 1, 4, 9 and 16, while at order // 6 , besides the simple poles at A 2 = 0, 1, 4, 9, 16 and 25, a double 
pole appears at N = 1, in agreement with the calculation of unitary integrals in reference ^3|. We can 

6 A "classical" Virasoro algebra, since there is no central charge. It is also known as a Witt algebra. 
7 We ignore the poles at unphysical negative values of N, which are completely irrelevant here. 
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check that our small- mass expansion of eqs. ( 5.13| ), ( |C.1| ) and ( |C.2| ) is correct by considering several 



special cases. For instance, if we take the mass matrix proportional to the identity M = m InxN 
all poles cancel, and, according to our normalization, we should obtain Z q = e mq . Its small-mass 
expansion is indeed in agreement with our results. This is one non-trivial check, for an arbitrary 
number of flavors and arbitrary value of q. We get another simple but non-trivial check for different 
masses by comparing our result with the small-mass expansion of the exact result [36] for the partition 



function Zq at fixed value of the number of flavors N. For example, for N = 2 with mass matrix 
M = diag(^i,/i2) and q = 0, the partition function is given by 

cr(2) ( \ sinhOi - /x 2 ) 

^qcdsCmi,/^) = — • (5.14) 

Mi — M2 

Expanding this expression for small masses, and comparing with eq. ( |5.13| ) and the corresponding 
expressions in Appendix C, we find perfect agreement. 

It is important to realize that all the poles at integer values of N cancel exactly once we specify 
whether N is even (q = 0) or odd (q = 1). Thus neither the even-iV nor the odd- A" partition functions 
have any physical poles at all. 



We have also checked that the expansion ( 5.13j ) is consistent with the expansion obtained from 



eq. ( |3.12| ) and the known small-mass expansion of the QCD4 partition function from ref. for 
Tr(M) = 0. It is instructive to follow the way the coefficients of the expansion of the QCD4 partition 



functions at u = ±1/2 combine to yield the corresponding coefficients for QCD3 as in eq. ( 5.13 ). One 



observes that the appearance of the v = ±1/2 factor on the right hand side of eq. ( 3.12| ) is related to 



the doubling in number of flavors on the left hand side. As a simple example, let us consider 

4(Nf ± v 




d N/) M = ^fl + 77^T + -l > (5-15) 



for Nj equal-mass flavors. This gives 



Z {Nf) (u)Z {Nf) (a) ~ 1 + ^u 2 I 1 + 1 1 + 
^u=l/2W^u=-X/2W ~ X + 2 ^ ^ 2^ + 1 + 2^-1 J + 



1 (2A r f) 2 o 
1+2(2^1^ + — (5 - 16) 



which coincides with our expansion ( 5.13 ) to this order (in agreement with equation ( |3.12| ) in the 



present normalization, where = 2Nf and q = 0.). 

For an even number of eigenvalues of the Dirac operator the effective finite-volume partition 
functions for QCD3, eqs. (|2.7|) and (|2,10 ), yield the same small- mass expansions as -2^qcd 3 (A^) 



(l/2)(Z ( q N) (MVZ) + Z^ q \MVT)) where q = 0, 1 for an even and odd number of flavors, respec- 
tively. Comparing the expansion (|5.13| ) with the small mass expansion of the QCD3 partition function, 
(ria=l det(.D + m a )) Yang-Mills ; we can read off the massless spectral sum rules for the QCD3 Dirac 
operator 
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N 2 — i 
~N{N r '- 

1-g 2 
N 2 -l 

(N 2 - 



1)' 



2 )(iV 2 -5q 2 -4) 



N(N 2 
(N 2 - 



l)(iV 2 -4)(iV 2 -9)' 



7 2 )(1 



iV(iV 2 -l)(iV 2 
g 2 [iV 4 (g 2 - 6) - 



-9)' 
8iV 2 (g 2 



5) +6(g 2 -4)] +3N 2 {N 2 -4) 



\ n ^ n m W , 



(N 2 



N 2 (N 2 -1)(N 2 - 
g 2 )(3iV 2 + 3g 2 -2iVV 



4)(iV 2 
12) 



9) 



N 2 (iV 2 -l)(iV 2 -4)(iV 2 -9) 




(iV 2 -g 2 )(24-10iV 2 + iV 4 



2„2\ 



6<f - A^g 



N 2 (N 2 - l)(N 2 -4)(iV 2 -9) 



(5.17) 



The sum rules obtained at order /j 6 are given in the Appendix B. It should be noticed that some sum 
rules given in ( 5.1 7| ) and (C.4) could not be derived had we assumed TrM 2k+1 = from the start. 
The point is that we only deduce sum rules for the massless theory. The quark masses only act as 
sources which should be as general as possible; otherwise some sum rules cannot be derived due to 
the vanishing of the corresponding source term. 

The first two massless sum rules agree with what was found in the original paper for even 
and odd flavors. We have also explicitly checked that the "diagonal" sum rules in ( |5.17D and ( |C.4j) 
agree with the results obtained from Random Matrix Theory by integrating the first two formulas of 
( [4.47 ). For the higher point functions a direct comparison with higher spectral correlation functions 
obtained from Random matrix theory becomes quite cumbersome. We have instead generated general 
relations among sum rules directly from the corresponding Random Matrix Theory via Schwinger- 
Dyson equations. For instance, from the identity 



,oo * i / i \ a * 

J n^£<% t A 2 ( {Aj)n 

J -°° 7 = 1 I Al \ k Ak ) p=l 



3 

we derive, after the rescaling Qk 

a+2 



2N f 



AfE 2 x2 
2 P 







E 



+ a 




(5.18) 



(5.19) 



0, 2, 4, as well as other similarly derived relations, are in 




We have checked that relation ( |5,19 ) for 
agreement with ( |5.17 ) and ( |C.4| ) . 

For even M, the spectral sum rules that are odd under Ck — * — Ck vanish identically as in the case 
of, e.g., (Sfc(l/Cfe)) = 0- This is a consequence of the symmetry Zq^ d ^(—M) = Z^} B ^(M). This 
symmetry is clear for ( 2.1C| ), and for (2.7) it follows from the following argument. In a small- mass 
expansion of ( 2.10| ) for each power of the matrix M we have the same power of the matrix Ts and the 
trace of odd powers of the mass matrix must be accompanied by the trace of odd powers of Ts which 
vanishes identically. 
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6 The Replica Limit of the QCD3 Finite Volume Partition Function 



In this section, starting from the replica limit of the small and large-mass expansions of the QCD3 
finite volume partition functions, we obtain the corresponding expansions of the partially quenched 
chiral condensate. 

6.1 Small-Mass Expansion 



Once we have the expansion of the QCD3 partition function, it should in principle be a simple matter 
to use the replica method to derive the partially quenched chiral condensate. However, in the case of 
the small-mass expansion this is not at all the case. The reason is the apparent poles in the expansion. 
As mentioned above, whether N is even or odd, all these poles cancel out exactly in the expansion 
of the partition function itself. But with the replica method we are forced to consider non-integer 
values of N = 2Nf + 2<n, where n is taken to zero. Then these so-called de Wit-'t Hooft poles [45] are 



significant: the analytic continuation cannot be performed across such poles. For larger and larger 
values of Nf the first pole is pushed to higher and higher order, but there will always be a finite order 
beyond which we cannot construct {/Uj}) from a small-mass expansion using the replica method. 

This issue was discussed at length in the case of QCD4 and we shall therefore be brief here. 

For simplicity, let us restrict ourselves to the case of massless physical fermions; the general case 
can be derived in a completely analogous manner. For a mass matrix M with 2Nf massless flavors 
and n replica flavors of mass +fi v paired with n replicas of mass — [i v we have N = 2n + 2Nf and, 

Ti(M k ) = 0, for odd k, 

Tr(M k ) = 2n/4 , for even k. (6.1) 
Using (|4.36j ), we formally find from the expansion of the partition function ( 5.13 ) and ( |C.2| ): 



EQy) m] - q 2 ) (AN 2 -q 2 ){hq 2 +4-4N 2 ) g 

£ 2N f {4N 2 - 1) ^ + 2N f (4N 2 - l)(4N 2 - 4)(4N 2 - 9) ^ 



+ 



2(4Nf - q 2 ) (4Nf - 4)(4Nf - 16) + 2lq A - Uq 2 4Nf + 140g 



+ • (6-2) 



2N f (4N 2 - 1)(42VJ - 4)(4iV| - 9)(47Vj - 16)(4iVj - 25) 



Here we have q = for integer Nf and q = 1 for half- integer Nf, and we always consider the case of 
an even number of Dirac eigenvalues. As stressed above, the above expansion is only meaningful up 
to the order at which a pole has to be crossed in order to take the replica limit. The result ( |6.2[ ) can 
now be compared with the expansion of the result obtained from Random Matrix Theory || 8 : 
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For positive replica mass. The condensate is odd in the mass, and appropriate absolute value signs have to be 
inserted for fx v < 0; for example ^ —* fx v \ix v \, etc.. 
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~^—m = 5) = g/^ - 105^ + 945^ + ' ' ' • (6 - 3) 

Comparing with the expression fl6.2|) obtained from the replica method we indeed find exact agreement 
precisely up to the order at which the replica method no longer can be applied, as a de Wit-'t Hooft 
pole has to be crossed. We nevertheless notice that, remarkably, all odd-order coefficients agree with 
those of Random Matrix Theory if one proceeds beyond the order at which the replica method ceases 
to be valid. This is perhaps less mysterious if one considers the fact that as Nf — * oo the replica 
expansion becomes exact to all orders (and this exact expansion has just odd powers). Because the 
replica method applied to the small-mass expansion only gives us the (correct) series up to the first 
de Wit-'t Hooft pole, we can only extract information about the microscopic spectral density of the 
Dirac operator up to that power in the mass which is free from such poles. 

6.2 Large-Mass Expansion 



In this subsection we use the replica method to calculate the large-mass expansion of the partially 
quenched chiral condensate. Because the expansion terminates, we are able to derive the exact mi- 
croscopic spectral density of the Dirac operator in QCD3 without relying on the assumption that it is 
the universal microscopic spectral density of Random Matrix Theory. We perform this calculation for 
one and two massless flavors by a direct replica calculation of the partition function (| 



We evaluate the partition function ( |2.7| ) with U an 2(Nf + n) x 2(Nf + n) unitary matrix. The 
mass matrix is taken to be diagonal with elements M^ + ^ f+n = —M^. Ultimately, we are interested 
in the case M\ = • • • = M^ f = 0, and the remaining replica masses equal to ±/x. The integral (2.7) 
can be rewritten as [ 25 1 , 



Z, 



(2N f +2n) 
QCD 3 



dX 



1 • • • 



dX Nf+n dU 1 dU 2 A 2 (A)e 



Tr M+ (Ui AUl +U 2 AUl ) 



(6.4) 



where U\ and U2 are (Nf + n) x (Nf + n) unitary matrices, M + is a diagonal matrix with diagonal 
elements given by M + kk = -Mfcfe- The diagonal matrix A contains the integration variables as diagonal 
entries. The Vandermonde determinant is denoted by A (A) = Y\k>i(Xk — Xi). The integrals over U\ 
and U2 are familiar Itzykson-Zuber integrals. The partition function is thus given by 



7 (2JV/+2n) 
"-QCD 3 



1 ^ et 2 e M fe Ai 

jx Nf+n A2(M) . 



(6.5) 



In the limit of Nf massless flavors and n quark masses equal to fi, the determinant can be rewritten 

as 



det e MkXl 



A(M) ~ ^U^kl^ll 



1 

Ai 



1 

X Nf-l 
A N f +n 



X Nf - 



\n— 1 _M-V/V>+n 
A N f +n e 



(6.6) 
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In general, this expression is rather complicated. However, the leading term in an expansion in 
can be obtained easily. We will first derive this term and then we will calculate the complete expansion 
for N f = 1 and N f = 2. 

There are two types of integrations: those that are accompanied by an exponential exp/iX^, and 
those that are not. The second type of integrals correspond to zero modes in a saddle point expansion 
and have to be evaluated exactly. The remaining integrals are evaluated perturbatively by writing 

dX k ■ ■ ■ = [ d\ k / dXk ■ ■ ■ ■ (6.7) 



— 1 J— OO J—OD 

and putting A& = 1 — Xf./fJ> in the integrals over the segment (— oo, 1], and = — 1 — x^/fj, in the 
integrals over the segment (— oo, —1]. In both cases the integral over the x k ranges over the segment 
[0, oo). Denoting by p the number of integrals with —1 as upper limit, the leading order contribution 
is obtained by keeping only the p = terms. The leading order term is then given by the diagonal 
term in the expansion of the square of the determinant. All terms are the same and we thus find 

Z QCD 3 +2n) = ,J+2N f n I , dXl ' ' ' I , dXN f ^ ( Al ' ' ' ' ' Xn f ) 



3 jJ n 2 +2N f n J x 

2a;jv f +1 2xjv f + n \ 2 



poo poo 

X 



dx N+l --- dx Nf+n e ZXN f +1 - ZXN f + "A 2 (x Nf+1 ,---,x Nf+n ). (6.8) 
Jo 

We recognize the partition function for the Laguerre ensemble which can easily be evaluated by means 
of the orthogonal polynomial method. However, we do not need this constant to find its leading order 
large-^i contribution to the resolvent. In the replica limit n — * it is simply given by 

4cT 2n ' P=0) = (6-9) 
so that the /i-dependence of the resolvent becomes 

S(/x) = 1 - ^ + 0(l/n 2 ). (6.10) 

Let us return again to the group integral fl2.7|) . For p / the integrand is invariant under a 
[U(n)/U(n — pi) x U{p\)\ x [U(n)/U(n — p-i) x U(p2)] submanifold of the coset. The integral over this 
manifold has to be performed exactly. Its volume is ~ n pl+P2 . In the replica limit of the resolvent 
only terms with p = p\ + p2 = 0, 1 contribute^. 

Next, we evaluate the integral ( |6.5| ) for Nf = 1. We expand this determinant with respect to its 
first row. Using the permutation symmetry of the integration variables the square of this expansion 
can be reduced to only two different contributions. The partition function thus simplifies to (up to an 
overall constant) 
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(n + l) e M(2A 2 +-+2A n+1 ) A 2 (A2) . . . ^ A;t+i) 



-n(n + i) e MAi+A 2 +2A 3 +-+2A„ +1 ) A(Al) Aa> . . . 5 An+1 )A(A 2 , A 3 , • • • , A n+1 ) 

(6.11 



9 It is important to stress that < pi , P2 < n and the reader may wonder how can we make sense of sums over all 
saddle points in the limit n — > 0. What we do in this section is equivalent to a presciption given in ^7|. This procedure 
has several caveats and is not mathematically rigorous J29| , but in all cases where it has been applied pfi|, p7| , |3l| ] it gives 
the correct result. 
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We will expand this partition function in powers of The first term in the integrand does not 

depend on Ai and therefore the integral over Ai of this term has to be performed exactly. It simply 
gives a factor 2. The remaining integrals are evolved perturbatively. We again split the integrals 
according to Q6.7| ) and put A& = ±1 — x k /fJ> depending on the upper boundary of the integral. For 
asymptotically large //, the integration range of the x k is thus given by [0, oo). In the first term the 
only nonvanishing possibility is pi = P2 = 0. The term with p\ = p2 = 1 is of order n 2 and does not 
contribute in the replica limit. However, in the second term we have the possibilities p\ = 1, p 2 = 
and pi = 0, P2 = 1 in addition the possibility p\ = p2 = 0. All other contributions vanish in the 
replica limit. We thus find that the asymptotic expansion of the QCD3 partition function is given by 
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(6.12) 



We evaluate these integrals by means of the orthogonal polynomial method ffij. The orthonormal 
polynomials corresponding to the exponential weight function e~ 2x are the Laguerre polynomials 

P k (x) = V2L k (2x). (6.13) 

By adding rows and columns the Vandermonde determinants can be expressed in terms of these 
orthogonal polynomials. Keeping only the leading order terms in 1/yu in the exponentially suppressed 
term, we arrive at the partition function 
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+ ■■■, (6.14) 



where the ellipses denotes higher order contribution in to the exponentially suppressed terms. 
Because the polynomials are normalized to one, the first integral simply gives n\. This can be easily 
seen by writing the determinants as a sum over permutations and performing the integrals by orthog- 
onality. In the second and third integral we have to treat the integration over x\ and X2 separately. In 
the same way as for the first integral, the remaining integrals just give (n — 1)1. To do the integrations 
over x\ and X2 we need the integral 

/•oo 

/ dxe- x P k (x) = (-l) k V2. (6.15) 
Jo 

In the third integral, only the terms that contain P n -i(x2) are nonvanishing. The integrals over x\ 
and X2 have to be performed separately, and, by orthogonality, the integrals over x 3 , ■ ■ • , x n+ i give a 
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factor (n — 1)!. We thus find 
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The valence quark mass dependence is thus given by 



E( M ) 



1-1+ 1 



-2ii 



li 2fi 2 2\x 2 ' 



(6.17) 



where we have calculated the replica limit of the exponentially suppressed term by writing l/(n — 1)! = 
n/n\. Notice that only the leading order term in the expansion in has been included in the 
exponentially suppressed term. 



The exact result for the microscopic spectral density for one massless flavor is given by 
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(6.18) 



The corresponding valence quark mass dependence is equal to 
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(6.19) 



in perfect agreement with the replica result. All higher order contributions to the exponentially 
suppressed term cancel! 



Let us now consider the case of two massless flavors. We again expand the determinant in ( |6 .6| ) 
and collect terms that result in the same value for the integral as we did in the case of one massless 
flavor. We find 
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(6.20) 



The integrations that are not accompanied by an exponential should be performed exactly, whereas 
the remaining integration variables are rescaled as A& = ±1 — x^j \x with sign determined by the upper 
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limit of the integration in ( |6,7| ). They are performed by expanding the determinants as sums over 
permutations. The integrations that are accompanied by the exponential e 2xk are easily evaluated by 
orthogonality, and give as result nl, (n — 1)! and (n — 2)!, respectively, for the three different types of 
terms in ( 6.2C| ). Summarizing, by keeping only the leading and first subleading 1/ fi corrections in the 
exponentially suppressed term, we find 
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In the last terms we have included a factor 2 to account for expanding A2 and A3 around +1 and — 1 
and vice versa. The last three terms represent the corrections to the Vandermonde determinant. 
They have been simplified by means of recursion relations for the Laguerre polynomials. The ellipses 
denote higher order corrections to the exponentially suppressed terms. To calculate the integrals we 
need the following result: 
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The partition function is thus given by 
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The sums over k and I are elementary, 
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In the replica limit this reduces to 
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which leads to the following u-dependence of the resolvent 
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We have not calculated the terms of order e 2/1 //i 4 that appear in the expansion of the partition 
function. By comparison with the exact result given by 
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(6.32) 



it turns out that such terms vanish in the replica limit (see also eq. (4.52)). To complete the calculation 
we will have to show that these terms, as well as higher order exponentially suppressed terms, vanish 
in the replica limit. The termination of the asymptotic series we have observed in this calculation is 
actually a more general phenomenon known as Duistermaat-Heckman localization [48, |2?J (sec [49] 
for a pedagogical discussion of this topic). 
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7 Conclusions 



In the domain where the Compton wavelength of the Goldstone modes is much larger than the size of 
the box, the QCD3 finite volume partition function for 2Nj flavors is given by a U (2Nf)/U(Nf) x U (Nf) 
unitary matrix integral. The main objective of this article has been to relate the spectrum of the QCD3 
Dirac operator to this partition function, and extensions of it with more flavors. An alternative starting 
point would have been the low-energy limit of the graded (or supersymmetric) generating function for 
the Dirac spectrum. This partition function is obtained by complementing the QCD partition function 
with both fermionic and bosonic "ghost quarks". It is invariant under an internal supersymmetry. 
In this way correlation functions of the QCD Dirac eigenvalues can be obtained in a mathematically 
rigorous way. However, we have here been interested in the question whether the usual low-energy 
effective partition function, without bosonic ghost quarks, already contains sufficient information to 
completely constrain the low-energy Dirac spectrum. 

This question has been answered affirmatively earlier within the context of Random Matrix Theory. 
The spectral fc-point correlation functions are given by the ratio of the usual partition function and a 
partition function with 2k additional flavors. These additional flavors can be related to the fermionic 
and bosonic ghost quarks of the supersymmetric generating function. A priori there is no reason why 
the spectral correlation functions of eigenvalues in Random Matrix Theory should coincide with those 
of the Dirac operator in the low-energy limit of QCD. Can we reach the same conclusion without 
relying on Random Matrix Theory? 

We have studied this question in three different ways. First, we have found exact relations between 
the QCD3 partition function and the QCD4 partition function analytically continued to half-integer 
topological charge. We have next added n replica quarks to the two theories and taken the replica 
limit. This has allowed us to derive explicit relations between the microscopic spectral density of the 
Dirac operator in QCD4 and the corresponding microscopic spectral density of the Dirac operator 
in QCD3 . These relations suggest the existence of not yet well-understood connections between the 
physics of QCD3 and QCD4, and are thus of independent interest. Second, the finite volume partition 
function of QCD3 satisfies Virasoro constraints. They can be solved recursively, providing us with 
the small mass expansion of the finite volume partition function, and a series of spectral sum rules. 
Third, we have directly analyzed the replica limit of the Unitary Matrix integral corresponding to the 
finite- volume partition functions of QCD3. In this way we have been able to derive the first few terms 
of the large mass expansion. Because the series truncate, this expansion turns out to be exact. 

A more general question we have addressed is the applicability of the replica method to the 
calculation of quenched averages. It was suggested recently that exact non-perturbative results can 
be obtained using an appropriate scheme of replica symmetry breaking. The general relations we 
have derived between the microscopic spectral densities of QCD3 and QCD4, though, follow from a 
completely straightforward replica limit. As an example, either by using previously obtained results 
for the QCD4 partition function or by the direct replica calculation of this paper we have reproduced 
exact results from the large mass expansion of the partition function. Possible caveats of the replica 
limit do not show up in our calculation. We do not know the reasons behind this observation, and 
it would be interesting to obtain a better understanding of the conditions for the applicability of the 
replica method. 

Generally, it has been expected that the replica method will not give exact results for spectral 
correlation functions but, at best, the coefficients of an asymptotic expansion. However, there are 
important cases where the series expansion turns out to truncate. We have witnessed a dramatic 
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example of this phenomenon in this article. However, there are still subtleties. The finite volume 
QCD4 partition function is well-defined for all integer values v of the topological charge. What is its 
analytical continuation in the complex z^-plane? A naive analytical continuations fails because of the 
presence of a shi7r^-term. The correct analytical continuation is obtained by changing the integration 
contour such that translational invariance is restored. We have not calculated for general Nf the 
order at which the large mass expansion terminates. We believe that it is important to get a better 
understanding of this question and its relation with Duistermaat-Heckman localization, but this is 
beyond the scope of this work. 

The replica method also yields the correct power-series expansion for small masses. As in QCD4 
the small mass expansion has de Wit-'t Hooft poles that limit the number of terms that can be 
obtained this way, but up to the order the replica method can be used all results agree exactly with 
earlier results from both Random Matrix Theory and the supersymmetric formulation. These three 
very different methods are thus obviously confirming each other. 

In conclusion, we have shown the existence of powerful relations between the QCD3 and the QCD4 
partition functions. The replica limit of these results provides us with exact non-perturbative relations 
between the spectral correlation functions of the Dirac operator in the two theories. It would be most 
interesting to derive similar relations from the supersymmetric formulation as well. 
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A Relations between QCD3 and QCD4 determinants 



In this Appendix we will prove the two eqs. ( [3.34 ) and ( [3.35 ) from two different lemmas, which then 
completes the proof of Theorem III eq. ( 3.27 ). 



Lemma 1: Let the matrices D({/ij}, £, lo), C({fii},C,uj) and S({/ij}) be defined as in eqs. (3.31) and 
(3.14), respectively. Then the following relation holds: 

(w + C) det D({ Mi }, C, w) + (w - C) det D({^}, -£ w) = (-1) JV /+ 1 2 7V /+ 1 det C({^}, C, w) det S({^}). 

(A.l) 

Proof: We will show that the left hand side factorizes into the right hand side by adding and 
subtracting rows and columns and expanding the determinants appropriately. 

Due to the block structure in the definition of the matrix D in eq. ( 3.31| ) we can add the upper 
blocks C and S to the lower ones and obtain after taking out factors of two, 
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When evaluating the determinant D({//j}, — £, u;) we can use that in the row depending on £ the 
functions C J_1 cosh^ -1 ^) are even under Q — ► — £ whereas the functions C"' -1 sinh^~ ^(C) are odd. In 
order to add the two determinants in eq. ( |A.l| ) we multiply the factors (u ± £) into the last column 
of the respective determinants. Next, we use that determinants differing by only one single column or 
row can be added. We use this to split the left hand side of eq. ( |A1| ) in the following way, where we 
explicitly display the last column, 



+ C)detD({// i },C,w) + (w-C)detD({/Zi},-C,w) 



2 N f 



det 



( C({ Mi }) 
C({Mi}) 



-s({M) 





C-?- 1 cosh( J - 1 )(C) -C i_1 smh^'- 1 )(C) 



cosh t''- 1 ) (a;) 



-a^sum^M 



+ det 



/ C({^}) -s({^» 
C({M) 
C^" 1 cosh^" 1 ) (C) -C j_1 sinh^'- 1 ) (C) 



3-1 



cos 



h a-D ( 



+ det 



+ det 



/ C({M) 
C({M) 

C J - 1 cosh( J - 1 )(C) 

/ C({M) 
C({Mi}) 

C J - 1 cosh^- 1 )(C) 



V 



COS 



-s({M) 



-aP'-^inh^- 1 ^^) 



-S({M) 


C^sinh^'-^C) 



(A.3) 

-^sinh^/)^) \ 


-C^ +1 sinh^)(C) 
-a; Ar /+ 1 sinh( Ar /)(a;) / 

-O^sinhW^) \ 


-wC^sinhWfC) 
-Cw^sinliW(w) / 



N r 



— w/^ sin 


sinh^C) 



(u N f sinh^^w) / 



In the next step we can add the first and the third determinant on the right hand side, since they only 
differ by the last but one row. And we can also add the second and forth determinant for the same 
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reason, after multiplying in the fourth determinant the last column and last but one row with (— 1). 
We find 



(w + C)detD(Oii},C,w) + (w-C)detD({^},-C,w) = 

( C({M) -s(M) 

C({M) 



2 jv>+i 



(A.4) 
\ 



det 



+ 



det 



+ 



det 







C' -1 coshG- 1 ) (C) -QN f +i sinhW) (£) 

V (J- 1 cosh^'" 1 ) (u) -<J~ X sinh^'- 1 ) (u) -w iV / +1 sinh W (w) / 

/ C({^}) -S({ M i}) -^f'sinhW(^) \ 

C({/ii}) 

C J '- 1 cosh^'- 1 )(C) 

V u^cosh^V) -uP'- 1 sinlr^ 1 V) / 

/ C({ Mi }) -S({M) -C^smhW)^) \ 

C({yUi}) 

-C^-^mh^'- 1 )^) 



where we have again split the determinant, resulting from the first addition, into two parts with respect 
to the last column, and taken out a factor of 2. The first determinant in the last equation is already 
part of the desired result. In the last column we can write sinhW =cosh( 7V / +1 ). If we then permute 
the last column Nf times to the left the first Nf + 2 columns only contain derivatives of cosh. We 
Laplace-expand the resulting determinant into products of (Nf + 2) x (Nf + 2) blocks and Nf x Nf 
blocks. After taking out all signs we obtain 



(-l)^ +1 2*/ +1 det 



/ C({M) 





■3-1 



cos 







\ (J- 1 coBk^iu) u; N f +1 coshW +1 )(cu) 



cosh^/ +1 )(C) detS({^}) . (A.5) 



What remains to be shown is that the two remaining determinants in eq. ( A.4 ) give rise to the same 



term as in eq. ( A.5 ), except that the last column vector is reading (fif f+1 cosh^ N f +1 \fii), 0, 0). In 



that way they will add up to the right hand side of Lemma 1, eq. ( A.l ) 



To proceed, we find that in the third determinant of eq. (A.4), the last term in the last column can 
be put equal to zero. To see this we expand the determinant with respect to the last column and note 
that the determinant multiplying —(^io Nf smlS N f \ui) vanishes, which can be easily seen by expanding 
it into products of (Nf + 1) x (Nf + 1) blocks and Nf x Nf blocks. Next, in both the second and third 
determinant in eq. (|A.4[) we take the common factors £ and u out of the last column, respectively, 
and multiply them into the last but one row. The two determinants are then equal up to precisely the 
last but one row and can thus be added. We arrive at 



2 N f +1 det 



( c(M) 



cos 



-s({M) 





cos 



htf- 1 )^) -^sinh^V) 






o J 



(A.6) 



It is straightforward to see that, when expanding with respect to the last column, and then into 
products of (Nf + 1) x (Nf + 1) blocks of C and Nf x Nf of S, the u multiplying the left side of the 
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last but one row can be multiplied into the right side of the last row 
2 N f +1 det 



( c«M) 

C({/*}) 







■j-i 



cos 



h o--i) (c) 





-^'sinh^'-^CC) 
-^'sinh^ -1 )(w) 



-^sinh(^)(/xO \ 






(A.7) 



As we have said before we claim that this is equal to the following product: 



(_l)JV/+i 2 JV /+ i det 



( c(M) 



/if /+1 cosh( JV / +1 )(/i J ) 



c 



■3-1 



COS 



3-1 



COS 



3h^- 1 )(C) 

h^ _1 )(w) 








detS({M) 



(Ai 



which will add up to eq. (A. 5) in order to give the desired right hand side of eq. (A.l). To see this we 
still have to do two more Laplace expansions. First of all we eliminate in eq. ( A.7 ) the upper block 
C by subtracting the lower one. Then, we permute the first column Nf times to the right to obtain 



2 N f +1 det 



( S(0*}) fj? f wihW(fM) \ 

^cosh^(/ij) cosh(/x.j) 
C J 'cosh^(C) cosh(C) C^sinh^^) 



\wicoBh.V\u) cosh(w) a^'sinh^" 1 ^^) 



(A.9) 



/ 



after taking out all signs. Our strategy will now be as follows. We will expand this determinant into 
products of Nf x Nf blocks and (Nf + 2) x (Nf + 2) blocks. The first blocks, after taking out common 
factors and rewriting coshV' = sinhv" - ^ will give rise to the desired prefactor det S(({/Ltj}). In doing 
so we obtain for eq. ( A.9 ) 



( _ 1) 7V /+ l 2 7V /+ l J det ^ cosh (i) (/U . 



Nf 

-£(-*) 

k=l 



det 



det 



' fi{ cosh^(/i^ 
. . . no k . . . 

y uji cosh^ (w) 

^ nl cosh^(^j) 

... no k . . . 
\ C J 'coBh^(C) 



S({M) 
del i cosh(C) C^'sinh^ 1 ^^) 
cosh(w) lu 3 ; sinh^' _1 ^(a;) 

cosh(^ fc ) 

del 1 S({A*i}) 

cosh(C) C sinh ' ; iO 

cosh(^) 

dd i S({/ii}) 

cosh(uj) uP ; sinh^' _1 ^(a;) 



/ , 








/Uj 7 sm 




(A.10) 





/^sinhW^ 




> . 



The first product gives already part of eq. ( |A.S| ) whereas the remaining contributions originate from 
the other terms. To see that we take out Yl i fii out of the first determinant and multiply the ith row 
of the second determinant with \ii (for i = 1, • • • , Nf). Then, we write cosh^ = sinh" -1 -' in the first 
and sinh^ -1 ' = cosh^) in the second factor to obtain 



(_)JV/+i 2 7V /+ 1 det S({/Xi}) det 



/ ^coshW(^) M f /+1 cosh( 7V / +1 )( W ) 

cosh(C) ^cosh^'^C) 
\ cosh(u;) i^'cosh^(w) 



(A.ll) 



We are still left with the sum over k in eq. ( A.lOj ). To simplify it further we take out common factors 
Hi^k from the first determinant as well as rewrite cosh^ = sinh^ -1 '. We then expand the second 
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determinant into products of 2 x 2 blocks and Nf x Nf blocks, where the 2x2 blocks consist of the 
first and last column 



N f 



2 N f+ l £ 
k,l=l 



k+l-1 



\i\ sinh^ x > (/ij 



det 



det 



... no k . . . 

I ^sinh^'- 1 )^) 
. . . no k . . . 

v C^smh^'-^C) 



cosh(/ife) fi^ } sinh.^ Nf ^ (^li) det 
cosh(^fc) ^f f smh.^ Nf \^i) det I 



' S({^}) 

v C^'sinh^- 1 )^) 

S({avi» 



2 N f +l 



N f 



£ (-l)^" 1 ^ I] A*i cosh(^) ^ +1 cosh^ +1 )( W ) 



fcj=l 



(detS({yUj^ fc },u;)detS({/ij^},C) - det S({/^ fe }, C) det S({^}, w)) 



2^+1 



Nf 



£(-l)*+'-iCw IJ Mi cosh(M fc ) Mr /+1 coBh( JV > +1 )( Ml )detS({A* i })detS({/z ¥ , )fc },C, 



AT, 



- II ^ cosh (^) Mf /+1 cosh^/ +1 )( W )detS({ W })detS({^M},C ; 

2 Ar /+ 1 detS({^}) 



x^(-i)^-vn M et 

(_l)^/+i 2 Ar /+ 1 detS({// i }) det 



cos 

COS 



h( W ) M f /+1 cosh^ +1 )( Mz ) 



detS({ ; u^ A:i /},C,w) 



/ COsh(^) /ifcosh^) M f^ +1 cOSh( Ar /+ 1 )(/X i ) 





V o 



C J 'cosh^')(C) 



0J° COS 



. (A.12) 



In the first step we have only taken out common factors £ and u> in order to be able to use the short 
notation S. In the second step we have used a property for antisymmetric products of determinants. 
A proof of the relations used can be found for example in Lemma 5 of [^5|. Moreover, we have to 
distinguish the cases k > I and k < I, whereas the terms with k = I drop out. In step three we 
have rewritten everything in terms of determinants. This is just to see that the double sum gives the 
Laplace expansion of the matrix in the last line. Here, we have again used that sinh^ -1 ) = cosh^'^ 
and multiplied the factor (oj Il j^fc.i M j m to the last determinant det S({ni^k,i}, C; w )- To summarize 
efforts, the last line of eq. ( [A. 12 ) together with eq. ( A.ll| ) sums up to eq. ( |A.S| ), as we wished to 



our 



show. So we see that starting with the left hand side of eq. (| 
which add up to the right hand side of the lemma eq. ( A.l[ ) 



we arrive at eqs. (|A.5|) and (A 



Lemma 2: With the same definitions as in Lemma 1 the following relation holds: 

{uj + C) det D({^.}, C, u)-(u- C) det D({^.}, w) = (-l) N f2 N f +1 det S({ M ;}, £ det C({^}). 

(A.13) 

Proof: The proof will go very much along the same lines as in Lemma 1. However, we will to arrive 
at a determinant of S of size (Nf + 2) instead of size Nf, which requires a different addition and 
subtraction scheme of rows and columns than in the proof of Lemma 1. 
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We start again from the definition of the matrix D in eq. ( |3.31 ) and this time subtract the upper 
blocks C and S from the lower ones, 



detD({/ii},C,w) = 2 JV /det 



/ o 

c({M) 



-s«M) 

S({M) 



\ 



C^'" 1 cosh( J '~ 1 )(C) -C^sinh^^C) 



(A.14) 



If we multiply the factors (u =t £) into the (Nf + l)-st column of the determinants and split them up 
the same way as before we obtain 



(w + C)detD({/ii},C, 



2 N f 



det 



w-C)detD({^},-C,w) = (A.15) 

( -8(M) \ 

C({M) c^f'coshW^) S({ W }) 

C^cosh^CO C^coshWjc) -C^sinh^-^C) 



+ det 



V c^" 1 cosh^'- 1 ) (w) w JV / +1 cosh W (u) 

( 

C({M) 



C/if'coshW( Mi ) 



-^-^inh^- 1 )^) / 

-s({M) 

S({M) 



C J '- 1 cosh(- ? - 1 )(C) wC^coshWCC) -^"^inh^- 1 )^) 



j-1 



cos 



h 0-D( 



det 



C({^}) 



Cw^/coshW^) 




V" 1 sinh^'- 1 V) / 



N 



S 



UJfl- COS 



-S({Mi}) 
S({M) 



C J '- 1 cosh(- ? - 1 )(C) -C^ +1 cosh^)(C) C i-1 sinh^'-^CC) 



cos 



det 



C(M) 



^ ) o; iV /+ 1 coshW (w) -u^'" 1 sinh^'" 1 ) (w) / 




-s({M) 

C^cosh^)^) S ({ Mi }) 



V 



C J '- 1 cosh (j - 1) (C) uC N t coshW(C) C^shm^ 1 ^) 
V u^'" 1 cosh^'- 1 ) (a;) -Qu N f cosh^ ) (w) 



-a^'-^inh^- 1 ^^ / 



In the next step we can again add the first and third determinant as well as the second and fourth 
determinant, after multiplying in the fourth determinant the last column and last but one row with 
(— 1). We obtain 



(u; + C)detD({ Mi },C,o;) 

= 2 N f +1 



det 



+ 



det 



C)detD({ w },-C,u/) = 

/ 

c({M) o 

C 7V/+1 cosh( 7V /)(C) 

V u> j ~ l cosh^'" 1 ) (u>) u N f +1 cosh W> (w) 

/ 

c({M) 



(A.16) 
\ 
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c({M) 





V 0J- 



cos 



3-1 



cos 



h^ _1 )(w) 



^coshW)( w ) 





CM^cosh^)^) 


Cc^coshW( 



-S({M) 
S({M) 
-C^sinhfr'-^C) 
-(J- 1 Baby- 1 ) (u) J 

-S({Mi}) \ 
S({M 

-C^sinh^^C) 
-S«Mi}) \ 

s({M) 







-u^sinh^V) / 
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where we have split the determinant resulting out of the first addition into two, with respect to the 
last column, and taken out a factor of 2. The first determinant in the last equation is already part of 
the desired result. In the (Nf + l)-st column we can write coshW =sinh( JV /+ 1 ). If we then permute it 
to the last column, the last Nf + 2 columns only contain derivatives of sinh. We then Laplace-expand 
the resulting determinant into products of Nf x Nf blocks and (Nf + 2) x (Nf + 2) blocks and obtain 

/ s({M) o \ 

(_)^/2 JV > +1 detC({^ i })det C^ 1 sinh^'- 1 ) (C) C Nf+1 sinh^+^C) (A. 17) 

\ ui-iembV-Viu) cu N f +1 sm^ N f +1 \uj) J 

after re-arranging all signs. What remains to be shown is that the two remaining determinants in 
eq. ( |A.16 ) give rise to the same determinant as in eq. ( |A.17 ), with the last column vector reading 
(H? f+1 sm^ N f +1 \fj,i),Q, 0) instead. 

We find again that in the third determinant in eq. ( |A.16| ), the last term in the (Nf + l)-st 
column can be set to zero following the same argument as before. Next, in both the second and 
third determinant in eq. ( A.16| ) we take the common factors £ and u> out of the (Nf + l)-st column, 
respectively, and multiply them into the last but one row. The two determinants are then equal up 
to precisely the last but one row and can thus be added, 

/ 



\ 



2 N f +1 det 



c({M) 

C J 'cosh( J - 1 )(C) 
(J- 1 coah^-^iu) 



cos 



S({/Xi}) 



(A.18) 



Along the lines of Lemma 1 it is again straightforward to see, that the factor uj multiplying the right 
side of the last but one row can be multiplied into the left side of the last row 







2 N f +1 det 



Nf 



COS 



-s({M) 

S({/*}) 



8111 



h^(0 



C({/ii}) /A 
C J 'cosh(^ 1 )(C) 
V o^cosh^^w) 

This remains to be shown to be equal to the following product: 

S(M) M f /+1 smhW +1 )(^) 



< 3 



(A.19) 



(_)^2 Ar / +1 detC({Mi})det 



3-1 
7-1 



Q 3 ' suit 



smr 



(A.20) 



ih^V) 

which then adds up to eq. ( [A. 17 ) to give eq. (A. 13). To proceed, we eliminate in eq. ( |A.19| ) the lower 
block S by adding the upper one. Next, we expand the determinant into products of (Nf + 2) x (Nf+2) 
blocks and Nf x Nf blocks. This time it will be the last blocks, that give us the prefactor det C(({ni}) 



-2 N f +1 I det 



C J 'cosh( J - 1 )(C) 
V t^cosh^ -1 )(w) 

/ 



Nf 



COS 



■If 
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C({M) 
V C^'sinh^- 1 )^) 

( ° 

C({M) 



Nf 







cos 



M f / cosh( Ar /)(/x 



sinh(C) det (^sinh^^)) 
sinh(w) / 

sinh(/j fc ) 


sinh(C) 
sinh(/j fc ) 


sinh(o;) 



(A.21) 
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det 



Hi sinh (j,) (/ij) 
. . . no k . . . 

uj 3 sinh^(a;) 

\i\ sinh^(/ij) 
... no k . . . 
C J 'sinhfr'>(C) 
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The first product gives part of eq. ( A.2C ) by taking Jli Mi ou t of the second determinant and multiplying 
the first Nf rows of the first determinant, upon rewriting cosh^ = sinh^ -1 ^, 



.)iv /2 iv /+ i det 



/ i4 sinh«) (jh) /v /+1 sinh( JV / +1 ) ( W ) 



smr 



h(C) v'sinlfiC) 
\ sinh(cj) tj- 7 sinh^^(tj) 



detC({M) . 



(A.22) 



Next, we tackle the sum over k in eq. (A. 21). Along the lines of Lemma 1 we expand the first 
determinant into products of Nf x Nf blocks and 2x2 blocks, 
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x (det C({/ij^ fc }, w) det C({/^/}, C) - det C({^ fc }, C) det C({/^}, w)) 
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^ sinh(^) M f /+1 smh( JV /+ 1 )( W )detC({ i u i })detC({/x iAfc },C ) u;) 
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-^(-l) fc+ 'C^Il ^ sinh (w0 Mr /+1 smh( Ar /+ 1 )(^)detC({// i })detC({ Mw },C,a;) 
fc</ j^M 

2 7V / +1 detC({^}) 



v 



xEf-^Il ^ ^ ( Sm Nf+1 . 

fc>i # fc ,j V smh( W ) fi t smW f + >(/*,) , 

sinh(/ij) ^sinh^^(^j) /+1 sinh^ +1 )(//j) 
(_l)^/ 2 JV /+ 1 detC({/i i }) det | C 5 'sinh^(C) 

i^'sinh^H 



h(^) ^ +1 sinh^+ 1 )( Atfc ) 



detC({^ M },C,w) 



(A.23) 



Here, we have performed the same steps as in the respective part of the proof of Lemma 1. Thus eqs. 
( |A.23j ) and ( A.22 ) sum up to eq. ( A.20| ), which together with eq. ( A.17| ) leads to the right hand side 
ofeq. (lAl3j ). 



B Relations between QCD3 and QCD4 from Random Matrix Theory 



In this subsection we show that the relations of Theorem III also can be derived from an entirely 
different point of view, based on the Random Matrix Theory representations of the QCD3 and QCD4 
finite- volume partition functions. 

Using the results of refs. f|, f§, we write the partition functions and (|2.12| ) as the large- AT 
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limits of the matrix integrals 



^QCD 3 



(W) 



2N f 



> Y[ det (0 + im/) exp -A/tr F(0 2 
/=i 



(B.l) 



for QCD3, and 



zl Nf \{ mi }) 



N f 



dW det + rrif) exp 
/=i 



-ytrV(^; 



(B.2) 



for QCD4. Here <j) is a Hermitian 2AA x 2AA matrix and 




(B.3) 



with TV an arbitrary complex matrix of size (A/" + ^) x A/". In both cases the Random Matrix Theory 
potentials V{cf) 2 ) are not constrained beyond giving convergent integrals with a non-vanishing spectral 
density at the origin, p(0) / ||. In eq. flB.l| ) the masses are taken to be pairwise grouped with 
opposite signs exactly as in the integral ([2.7j), The limit M — > 00 is taken in the two Random Matrix 
Theory partition functions so that \X{ = miirp(0)J\f and pn = miirp(0)2Af, respectively, are kept fixed. 
In this limit the partition functions Zqq^ and Zqq^ become proportional to each other, with a 

Similarly, also the partition functions zj) Nf ^ and Z^ f ^ 



/Uj-independent proportionality constant [ 
become proportional Below we derive relations between the unitary matrix integrals (2.7) and 
( [2.12 ) from their Random Matrix representation ( B.l ) and (|B.2j ), respectively. The key ingredient is 



to use the standard orthogonal polynomial technique of Random Matrix Theory, and next rewrite all 
expressions in terms of the partition functions themselves [40, g]. Additional details can be found in 
ref. m. 



We begin by expressing the integrals ( |B.1| ) and ( B.2[ ) in terms of eigenvalue representations. This 
permits us to consider non-integer values of v through analytical continuation. The starting point is 
a relation between the orthogonal polynomials of chUE (QCD4) and those of the UE (QCD3): 



I, chUE 



-1/2) 



({m i };z' 2 



R 



(2N f ) 



21, UE 



(B.4) 



This simple identity is easily derived from the definition of the two Random Matrix Theories involved 
Since in ref. [|39| no explicit use was made of the measure, eq. ( |B.4[ ) also holds in the massive 
case. On the right hand side the same Nf masses appear, but doubled into pairs with opposite sign. 

Instead of using directly the orthogonal polynomials P/({mj}; A), it is convenient to work with the 
"wave functions", 

(B.5) 



Vi({mi};\) = ^to({m i },X)Pi({m l };\) 



where w({m,j; A) is the measure function (so that the wave functions f;({mj}; A) are orthogonal with 
respect to a weight of unity). It follows from eq. (3.11) of ref. |3^] that one can use the wave functions 
of the UE inside the Christoffel-Darboux identity for the chUE kernel. We then find the following 
identity M]: 
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K 
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z 2 — u> 



2 V V 2AT, UEl z J ™ V 2AT-1, UE 
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>)*§£'ue(«o 



1 c. 



2 (z 



2Af 

w) 



\^2Af, VE\ z )^2Af-l, VE\ W > ~ W 2Af-l, XJE\ z )^2Af, VE\ W ) 



1 C, 



(2N f ) 
2Af 



2{z + w) 



2A/\ UE( z i V 2Ar-l, UeI^J + W 2A^-1, \]EK z )^2Af, UE^ 
(2JV» 



v(2JV/) 



S^UE^.™) + K 2N, f XJE(- z > w ) » 



(B.6) 



where for clarity we have not explicitly indicated the mass dependence of the wave functions. The 
coefficients c^fj are defined in ref. |3j|. We have also used that the parity of the polynomials Pi(X) 
for the UE is equal to (— 1) . The different kernels can all be expressed in terms the corresponding finite 
volume partition functions [{40| , ||, as we will see below. In fact, the relation (BJ3) simply expresses 
Lemma I of Appendix A, divided on each side by the result of Theorem I. 

In order to arrive at a more general result, we will also derive an equation involving Z^j\ / 2 ({/•*})• 



To do so we start with the kernel K 



(N f , u=+l/2) 
Af, chUE 



(z 2 ,w 2 ) with v = +1/2 in contrast to eq. (I 



Here, 



we will make use of the following relation among the wavefunctions of the UE: 



(2N f ,2a+2) ( , 
w 2«, UE \ z ) 



sign(z)^ {2 ^ ] E (z) 



(B.7) 



This relation was shown in |^] for any number of massless flavors 2a, but it can easily be generalized 
to include any additional number of 2Nf (massless or massive) flavors as well, since they are merely 
spectators. The sign function takes into account that the left hand side is an even function. Employing 
again the relation eq. (|B,4|) we obtain 
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sign(z)sign(w) 
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(B.8) 



In the first step we have used the relation ( |B.7| ) and in the second step the recursion relation XPi(X) 
q + iP; + i(A) + qP;_i(A). With the same rewriting as in eq. ([B.6D we arrive at 



(N f , u=+l/2) ( 2 2, 
Af, chUE \ z i w ) 



(2N f +2) (2N f +2) 
c 2Af c 2Af-l 



1 

2 
+ 



c 2Af 
1 

[z — w) 

1 

(z + to) 



sign(z)sign(u;) 



(B.9) 



(2N f ) 
2A/+1, UE 



~ w 2Af+l, UEV 



(2JV/) 
2A/\ UE 



(2W/) 
2A/+1, UE 



(2N f ) 
2Af, UE 



(2iV/) 
2A/+1, UE 



H. 
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(2JV/+2) (2JV/+2) 

= sign(z)sign(w) C2A/ " (2jV/) (#££+J, ue(*> «0 ~ ^SK, ue(-«> ™)) • 

c 2Af C 2AT+1 

In the large— jV limit the mass dependence of the recursion coefficients c^n^ becomes subleading and 
these coefficients simply cancel out. Again, this relation has a direct interpretation in terms of finite 
volume partition functions: it is Lemma II of Appendix B divided on each side by the relation of 
Theorem I. 



Adding and subtracting eqs. (R6) and (|B.9| ) we finally arrive at 

k^\m,±m = 1{k^T~ 1 ^{ J hUW) ± ^fflr +1/2) ({M,C 2 ,- 2 )) (b.io) 

in the microscopic large-AA limit, for £, uj > 0. One of these (the one with "+" signs) was previously 
known only in the massless case, where it can readily be derived from the generalized Laguerre ensemble 



We now make the following observation. From the kernel relation (|B.10| ) it follows that all spectral 
correlation functions of the two Random Matrix Theories ( |B.1| ) and (|B.2| ) are related in a fairly simple 
way. If one believes that these relations extend to the microscopic spectral correlators of the Dirac op- 
erators in QCD3 and QCD4, then one obtains all spectral correlation functions Pqcd 3 (Cii ■ ■ ■ > Cfei {fa}) 
in QCD3 from the kernel relevant for QCD4. The most obvious relation that follows from ( |B.10| ) is 
the one for the microscopic spectral densities, 

p®C;{M) = \ [pSC +1/2) (C;^» + P$gr 1/Si • (B.11) 

This relation is here derived within the context of Random Matrix Theory. As shown in this paper, 
it also follows directly from the field theory formulation, through the replica method. 

We note that based on the exact relationship of two-point correlation functions to the Random 
Matrix Theory partition functions |40| , ^] (which in turn are proportional to the finite-volume field 
theory partition functions), eq. (B.10) implies the following identity: 



This particular relation is simply the identity of Theorem III divided on each side by the identity 
of Theorem I, both of which were proven in section 4.1. The right hand side involves a sum or 
difference of two different ratios of partition functions that cannot easily be combined. The identity 
as it stands is therefore not suitable for deriving, by means of the replica method, identities among 
spectral correlators of the Dirac operators in the two theories. We note, however, that in the limit, 
£, uj — ► it is evidently consistent with the general identity ( |3.12 ). In this limit the second term on the 



right hand side of ( B.12| ) vanishes. Although the first term is also multiplied by ((w) 1//2 , the product 



of the two remains finite. Indeed, invoking flavor-topology duality [^g], we can express the product in 
terms of a partition function with two massless particle less, at the cost of increasing v, 



g£ g g) ({M,o,Q) * i_fgg 3/2 ({M) 



^ T N T • (B.13) 
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By applying flavor-topology duality, the relation ( |B.13| ) also follows from the more general identity 
( p. 12 ) given in Theorem I alone. To see this, we take the ratio of eq. ( 3.12j) for Nf + 1 and Ay and 



then send the extra mass to zero. The factors cancel after the flavor-topology shift ( |3.11|) 

on the partition functions with v = —1/2 and v = +1/2 and we arrive at eq. ( B.13 ). 



C Partition Function and Spectral Sum Rules to Sixth Order 



In this Appendix we give the terms of order /j 5 and // 6 in the small expansion of the QCD3 partition 
function Zq N \ The terms up to order fj, 4 were already given in eq. ( 5.13 ). For the terms of 0(fi 4 ) we 
find 



^q 111 5 



2q(N 2 - q 2 )(20 - 3N 2 + 7q 2 



-TrM 



5A(A 2 - 1)(A 2 - 4) (A 2 - 9) (A 2 - 16) 
q(N 2 - q 2 ) [N 4 + A 2 (4 - 5q 2 ) + 24(q 2 - 4)] * 

4A 2 (A 2 - 1)(A 2 - 4) (A 2 - 9)(N 2 - 16) 
q(N 2 - q 2 )(N 2 q 2 + 12q 2 - A 4 + 12A 2 - 48) 



(C.l) 



+ 
+ 



3iY2(AT2 _ 1)(AT2 _ 4 )(AT2 _ 9)^2 _ 16 ) 
q(N 2 - q 2 )(4 - q 2 ) 



TrM 3 TrM 2 



3A(A 2 - 1)(N 2 - 4)(iV 2 - 16) 
q{N 2 - q 2 ){N 4 - N 2 q 2 - 18 A 2 + 2q 2 + 



TrM 3 (TrM) 2 

TrM (TrM 2 ) 2 



8A(A 2 - 1)(A 2 - 4) (A 2 - 9) (A 2 - 16) 
q(N 2 - q 2 )(N 4 q 2 - UN 2 q 2 + 24q 2 - 3N 4 + 40iV 2 - 96) 

12N 2 (N 2 - 1){N 2 - 4) (A" 2 - 9)(N 2 - 16) 
q [N 4 (q 4 - 10q 2 + 15) - 20A 2 (g 4 - 9q 2 + 11) + 6(13<? 4 - 100g 2 + 96)] 



TrM 2 (TrM) ; 



120 A( A 2 - 1)(A~ 2 - 4) (A 2 - 9) (A 2 - 16) 



(TrM) £ 



At order fj 6 we find 



^q I A* 6 



AiTtM 6 ~ 1 (TiM) 1 + SfcTrM 6_2fc (TrM 



2\k 



k=l 



+ 



J2 Q (TrM^TrM 3- *) 2 + D (TrMTrM 2 TrM 3 ) . 



i=l 



where 
A = 
Ax = 
A 2 = 
A 3 = 



(N 2 - q 2 ) [{N 2 - 4) (AT 2 - 16) + 140g 2 - 14A 2 g 2 + 21 g 4 ] 

3A(A 2 - 1)(A 2 - 4)(A 2 - 9)(A 2 - 16)(A 2 - 25) ' 
2(A~ 2 - q 2 ) [-5(A 2 - 4) (A 2 - 16) + q 2 (-200 - 25A~ 2 + 3A 4 + 70g 2 - 7N 2 q 2 )] 

5N 2 (N 2 - 1)(N 2 - 4) (A 2 - 9){N 2 - 16) (A 2 - 25) 
(1 - q 2 )(A 2 - q 2 ) [{N 2 - 4) (A 2 + 23) + 13g 2 - 5A 2 q 2 ] 
8N(N 2 - 1) 2 (N 2 - 4) (A 2 - 9) (A 2 - 25) ' 

-(A 2 - q 2 ) 

18A 2 (A 2 - 1) 2 (A 2 - 4)(A 2 - 9)(A 2 - 16)(A 2 - 25) 



1541AV + 356A 4 g 2 - 15AY - 60g 4 + 229A 2 g 4 - 51A 4 g 4 + 2N & q 4 



(C.2) 



(A 2 - 4) (A 2 - 16) (9A 2 + 15) 
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^4 

A 5 



(N 2 - q 2 ) [3(N 2 - 4)(iV 2 - 16) + g 2 (-280 + 136iV 2 - 6iV 4 + 38g 2 - 24N 2 q 2 + N 4 q 2 )] 



48N(N 2 - 1) 2 (N 2 - 4) (TV 2 - 16) {N 2 - 25) 



-i 



15iV 2 (4 - N 2 )(N 2 - 9)(N' Z - 16) 



B 2 

Ci 
C 3 

D 



72QN 2 {N 2 - 1) 2 {N 2 - 4)(iV 2 - 9){N 2 - 16) (TV 2 - 25) 

+5760g 2 - 29416iV 2 g 2 + 13756iVV - 1545iV 6 g 2 + 45iVV - 2400g 4 + 14290iVV 
-58507V 4 q 4 + 575iVV - 15iVV + 240g 6 - 12587V 2 q 6 + 458iVV - 41iV 6 g 6 + N 8 q 6 

~( N2 - ^ [(40 - 17iV 2 + iV 4 )(iV 2 - 4)(iV 2 - 16) 

8N 2 {N 2 -1) 2 {N 2 -4){N 2 -9)(N 2 -16)(N 2 -25)[ { ^^+^ A iv V[ iy ioj 

+ g 2 (-800 - 484iV 2 + 907V 4 - 67V 6 + 40g 2 + 75N 2 q 2 + 5iV 4 g 2 ) 
(iV 2 - q 2 ) 



(N - 27 N z + 98) (TV 2 - 4) (TV - 16) 



48iV(iV 2 - l) 2 (iV 2 - 4)(N 2 - 9)(iV 2 - 16) (TV 2 - 25) 
+ <? 2 (-3480 - 164iV 2 + 467V 4 - 2iV 6 + 358<? 2 + N 2 q 2 + N 4 q 2 ) , 

(N 2 - q 2 ){q 2 - 1) [{N 4 - 15N 2 - W)(N 2 - 4) + g 2 (-10 + 3iV 2 - N 4 ) 



a 



32N 2 (N 2 - 1) 2 (N 2 - 4)(iV 2 - 9)(N 2 - 25) 
2(iV 2 - q 2 ) 



(15 - 3JV 2 )(iV 2 - 4)(iV 2 - 16) 



9iV 4 (iV 2 - l) 2 (iV 2 - 4)(iV 2 - 9)(iV 2 - 16) (TV 2 - 25) 

+ <? 2 (-600 + 3iV 2 - 47V 4 + N 6 + 90q 2 - 29N 2 q 2 - N 4 q 2 ) 

{N 2 - g 2 ) [3(N 2 - 4)(iV 2 - 16) - g 2 (180 - 32iV 2 + 2iV 4 - 13g 2 - 2N 2 q 2 ) 
6iV(iV 2 - 1) 2 (A^ 2 - A)(N 2 - 16) {N 2 - 25) 



(C.3) 



From ( |C.2| ) and ( |C.3| ) we derive the following additional spectral sum rules for an even number of 
Dirac eigenvalues: 



\^ Cn , 
\ n ^ n / m ^rn / 



6 Ac 



48 A A 



-6! A, 



X ~F~ X 72~ X ?3 

\ n ^ n m W p ^p i 



1 



1 



1 



1 



Cn 
1 



6 A 



8A 2 , 



1 



Sri 



1 



2\ 



( X 72 X/ 74" 

\ n m >m i 



48^2, 



= -18A 3 , 
-18iV 2 C 3 
> = -8 Si. 

We remind the reader that g = , 1 for even and odd number of flavors respectively. 



32 d, 



(C4) 
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